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Abstract 

We review and extend in several directions recent results on the "asymptotic safety" approach 
to quantum gravity. The central issue in this approach is the search of a Fixed Point having 
suitable properties, and the tool that is used is a type of Wilsonian renormalization group equation. 
We begin by discussing various cutoff schemes, i.e. ways of implementing the Wilsonian cutoff 
procedure. We compare the beta functions of the gravitational couplings obtained with different 
schemes, studying first the contribution of matter fields and then the so-called Einstein-Hilbert 
truncation, where only the cosmological constant and Newton's constant are retained. In this 
context we make connection with old results, in particular we reproduce the results of the epsilon 
expansion and the perturbative one loop divergences. We then apply the Renormalization Group 
to higher derivative gravity. In the case of a general action quadratic in curvature we recover, 
within certain approximations, the known asymptotic freedom of the four-derivative terms, while 
Newton's constant and of the cosmological constant have a nontrivial fixed point. In the case 
of actions that are polynomials in the scalar curvature of degree up to eight we find that the 
theory has a fixed point with three UV— attractive directions, so that the requirement of having a 
continuum limit constrains the couplings to lie in a three-dimensional subspace, whose equation is 
explicitly given. We emphasize throughout the difference between scheme-dependent and scheme- 
independent results, and provide several examples of the fact that only dimensionless couplings 
can have "universal" behavior. 



*Electronic address: a.codello@gmail.com 
^Electronic address: percacci@sissa.it 
* Electronic address: rahmedc@sissa.it 



2 



I. INTRODUCTION 



It is well known that general relativity can be treated as an effective quantum field theory 
. This means that it is possible to compute quantum effects due to graviton loops, 
as long as the momenta of the particles in the loops are cut off at some scale. For example, 
in this way it has been possible to unambiguously compute quantum corrections to the 



Newtonian potential {2]. The results are independent of the structure of any "ultraviolet 
completion" , and therefore constitute genuine low energy predictions of any quantum theory 
of gravity When one tries to push this effective field theory to energy scales comparable to 
the Planck scale, or beyond, well-known difficulties appear. It is convenient to distinguish 
two orders of problems. The first is that the strength of the gravitational coupling grows 
without bound. For a particle with energy p the effective strength of the gravitational 
coupling is measured by the dimensionless number \/Zx, with G = Gp 2 . This is because the 
gravitational couplings involve derivatives of the metric. The consequence of this is that if 
we let p — > 00, also G grows without bound. The second problem is the need of introducing 
new counterterms at each order of perturbation theory. Since each counterterm has to be 
fixed by an experiment, the ability of the theory to predict the outcome of experiments is 
severely limited. 

Could a Quantum Field Theory (QFT) get around these obstacles? Recall that Newton's 
constant, as any coupling constant in a QFT, must be subject to Renormalization Group 
(RG) flow. It is conceivable that when p — * 00, G(p) ~ p~ 2 , in which case G would cease to 
grow and would reach a finite limit, thereby avoiding the first problem. If this is the case, 
we say that Newton's constant has an UV Fixed Point (FP). More generally, if we allow 
the action to contain several couplings gi with canonical mass dimension di, we say that the 
theory has a FP if all the dimensionless parameters 

9i = Qik^ (1) 

tend to finite values in the UV limit 1 . This particular RG behavior would therefore solve 
the first of the two problems mentioned above and would guarantee that the theory has a 
sensible UV limit. 



1 Strictly speaking only the essential couplings, i.e. those that cannot be eliminated by field redefinitions, 
need to reach a FP. See [s| for a related discussion in a gravitational context. 
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In order to address the second problem we have to investigate the set of RG trajectories 
that have this good behavior. We want to use the condition of having a good UV limit as 
a criterion for selecting a QFT of gravity. If all trajectories were attracted to the FP in 
the UV limit, we would encounter a variant of the second problem: the initial conditions 
for the RG flow would be arbitrary, so determining the RG trajectory of the real world 
would require in principle an infinite number of experiments and the theory would lose 
predictivity. At the other extreme, the theory would have maximal predictive power if there 
was a single trajectory ending at the FP in the UV. However, this may be too much to ask. 
An acceptable intermediate situation occurs when the trajectories ending at the FP in the 
UV are parametrized by a finite number of parameters. A theory with these properties was 
called "asymptotically safe" in Q]. 

To better understand this property, imagine, in the spirit of effective field theories, a 
general QFT with all possible terms in the action which are allowed by the symmetries. 
We can parametrize the (generally infinite dimensional) "space of all theories", Q, by the 
dimensionless couplings gi. We assume that redundancies in the description of physics 
due to the freedom to perform field redefinitions have been eliminated, i.e. all couplings 
are "essential" (such couplings can be defined e.g. in terms of cross sections in scattering 
experiments). We then consider the Renormalization Group (RG) flow in this space; it is 
given by the beta functions 

If there is a FP, i.e. a point with coordinates g^ such that all fli(g*) = 0, we call C its 
"critical surface" , defined as the locus of points that are attracted towards the FP when 
k — > oo 2 . One can determine the tangent space to the critical surface at the FP by studying 
the linearized flow 

k ^ 1 ^ =M lJ (~g J -~9i*), (3) 



where 

d/3i 



M, 



dg 



(4) 



2 RG transformations lead towards lower energies, and the trajectories lying in C are repelled by the FP 
under these transformations. For this reason, C is also called the "unstable manifold" . Since we are 
interested in studying the UV limit, it is more convenient to study the flow for increasing k. 
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The attractivity properties of a FP are determined by the signs of the critical exponents 
defined to be minus the eigenvalues of M. The couplings corresponding to negative 
eigenvalues (positive critical exponent) are called relevant and parametrize the UV critical 
surface; they are attracted towards the FP for k — > oo and can have arbitrary values. 
The ones that correspond to positive eigenvalues (negative critical exponents) are called 
irrelevant; they are repelled by the FP and must be set to zero. 

A free theory (zero couplings) has vanishing beta functions, so the origin in Q is a FP, 
called the GauBian FP. In the neighborhood of the Gaufiian FP one can apply perturbation 
theory, and one can show that the critical exponents are then equal to the canonical dimen- 
sions ($j = di), so the relevant couplings are the ones that are power-counting renormalizable 
3 . In a local theory they are usually finite in number. Thus, a QFT is perturbatively renor- 
malizable and asymptotically free if and only if the critical surface of the GauBian FP is 
finite dimensional. Points outside C flow to infinity, or to other FP's. A theory with these 
properties makes sense to arbitrarily high energies, because the couplings do not diverge in 
the UV, and is predictive, because all but a finite number of parameters are fixed by the 
condition of lying on C. 

If the GauBian FP is replaced by a more general, nontrivial FP, we are led to a form of 
nonperturbative renormalizability. It is this type of behavior that was called "asymptotic 
safety". It is completely equivalent, in the language of critical phenomena, to the assumption 
that the world is described by a "renormalized trajectory" emanating from some generic FP 
{15I ] . In general, studying the properties of such theories requires the use of nonperturbative 
tools. If the nontrivial FP is sufficiently close to the GauBian one, its properties can also 
be studied in perturbation theory, but unlike in asymptotically free theories, the results of 
perturbation theory do not become better and better at higher energies. 

In order to e S tab ll8 h whether gravity has this type of behavior, several authors fl fl 
I2J applied the e expansion around two dimensions, which is the critical dimension where 
Newton's constant is dimensionless. The beta function of Newton's constant then has the 
form 

/3 d = eG+ B X G 2 , (5) 
where G = Gk~ e and B 1 < , so there is a FP at G — —e/B\ > 0. Unfortunately this 



The behavior of dimensionless or marginal couplings may require a more sophisticated analysis. 
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result is only reliable for small e and it is not clear whether it will extend to four dimensions. 

In the effective field theory approach (in d — 4), Bjerrum-Bohr, Donoghue and Holstein 
have proposed interpreting a class of one loop diagrams as giving the scale dependence of 
Newton's constant [4|. They calculate 



G(r) = Go 



167 G 



30tt r 2 

where r is the distance between two gravitating point particles. If we identify k = 1/ar, 
with a a constant of order one, this would correspond to a beta function 

Pa = 2G - a 2 ^G 2 . (6) 
G 15tt v ' 

This beta function has the same form as ([5]) in four dimensions, and, most important, 
the second term is again negative. This means that the dimensionful Newton constant G 
decreases towards lower distances or higher energies, i. e. gravity is antiscreening. This is the 
behavior that is necessary for a FP to exist, and indeed this beta function predicts a FP for 

= i§^2- This calculation was based on perturbative methods and since the FP occurs at 
a not very small value of G, it is not clear that one can trust the result. What we can say 
with confidence is that the onset of the running of G has the right sign. Clearly in order to 
make progress on this issue we need different tools. 

In this paper we will discuss the application of Wilsonian renormalization group methods 
to the UV behavior of gravity. In section II we will introduce a particularly convenient tool, 
called the "Exact Renormalization Group Equation" (ERGE) which can be used to calculate 
the "beta functional" of a QFT. Renormalizability is not necessary and the theory may have 
infinitely many couplings. In section III we illustrate the use of the ERGE by calculating 
the contribution of minimally coupled matter fields to the gravitational beta functions. In 
this simple setting, we will review the techniques that are used to extract from the beta 
functional the beta functions of individual couplings, emphasizing those results that are 
"scheme independent" in the sense that they are the same irrespective of technical details 
of the calculation. In section IV we apply the same techniques to the calculation of the 
beta functions for the cosmological constant and Newton's constant in Einstein's theory in 
arbitrary dimensions, extending in various ways the results of earlier studies 
We also show that the FP that is found in four-dimensional gravity is indeed the continuation 
for e — > 2 of the FP that is found in the 2 + e expansion. We compare various ways of defining 
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the Wilsonian cutoff and find the results to be qualitatively stable. In sections V and VI 
we make connection with old results from perturbation theory. In section V we rederive 
the 't Hooft-Veltman one loop divergence from the ERGE and we show it to be scheme- 
independent. We also discuss why the Goroff-Sagnotti two loop divergence cannot be seen 
with this method and we discuss the significance of this fact. In sections VI and VII we 
consider higher derivative gravity. In section VI we derive the existence of the FP in the most 
general truncation involving four derivatives at one loop, and we highlight the differences 
between the Wilsonian procedure 3] and earlier calculations. In section VII we consider 
higher powers of curvature, restricting ourselves to polynomials in the scalar curvature. We 
give more details of our recent calculations [ljj] and extend them to polynomials of order 
eight. In section VIII we assess the present status of this approach to quantum gravity and 
discuss various open problems. 

II. THE ERGE AND ITS APPROXIMATIONS 

The central lesson of Wilson's analysis of QFT is that the "effective" (as in "effective field 
theory" ) action describing physical phenomena at a momentum scale k can be thought of as 
the result of having integrated out all fluctuations of the field with momenta larger than k 
[15I ] . At this general level of discussion, it is not necessary to specify the physical meaning of 
k: for each application of the theory one will have to identify the physically relevant variable 
acting as k 4 . Since k can be regarded as the lower limit of some functional integration, we 
will usually refer to it as the infrared cutoff. The dependence of the "effective" action on k 
is the Wilsonian RG flow. 

There are several ways of implementing this idea in practice, resulting in several forms of 
the RG equation. In the specific implementation that we shall use, instead of introducing 
a sharp cutoff in the functional integral, we suppress the contribution of the field modes 
with momenta lower than k. This is obtained by modifying the low momentum end of the 
propagator, and leaving all the interactions unaffected. We describe here this procedure for 
a scalar field. We start from a bare action S[(j)], and we add to it a suppression term ASk[4>] 
that is quadratic in the field. In flat space this term can be written simply in momentum 



4 In scattering experiments k is usually identified with some external momentum. See 16| for a discussion 
of this choice in concrete applications to gravity. 
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space. In order to have a procedure that works in an arbitrary curved spacetime we choose 
a suitable differential operator O whose eigenfunctions <p n , defined by (D<p n = A n <£>„, can be 
taken as a basis in the functional space we integrate over: 

n 

where <fi n are generalized Fourier components of the field. (We will use a notation that is 
suitable for an operator with a discrete spectrum.) Then, the additional term can be written 
in either of the following forms: 

AS k [<P] = l -j dx <j>(x) R k (0) <j>(x) = ^0^ fe (A„). (7) 

The kernel Rk{0) will also be called "the cutoff". It is arbitrary, except for the general 
requirements that R k {z) should be a monotonically decreasing function both in z and k, 
that R k (z) — > for z ^> k and R k (z) 7^ for z <C k. These conditions are enough to 
guarantee that the contribution to the functional integral of field modes <p n corresponding 
to eigenvalues A n <C k 2 are suppressed, while the contribution of field modes corresponding 
to eigenvalues A n ^> k 2 are unaffected. We will further fix R k (z) —> k 2 for k — *■ 0. We define 
a /c-dependent generating functional of connected Green functions by 

e~ Wk[J] = J Dcf>exp {-S[(f)] - AS k [0] - j ' dxJ<j> 

and a modified /c-dependent Legendre transform 

T k [<t>} = W k [J] - j dxJ<P- AS k [<P] , 

where ASk[4>] has been subtracted. The functional T k is sometimes called the "effective 
average action", because it is closely related to the effective action for fields that have 
been averaged over volumes of order k~ d (d being the dimension of spacetime) 171 ] . The 
"classical fields" 5W k /5J are denoted again for notational simplicity. In the limit k — > 
this functional tends to the usual effective action the generating functional of one- 

particle irreducible Green functions. It is similar in spirit to the Wilsonian effective action, 
but differs from it in the details of the implementation. 

The average effective action r fe [</>], used at tree level, gives an accurate description of 
processes occurring at momentum scales of order k. In the spirit of effective field theories, 
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we shall assume that Tk exists and is quasi-local in the sense that it admits a derivative 
expansion of the form 

oo 

r k (<P,9 l ) = J2J2^\k)Oi n) (<P) , (8) 

n=0 i 

where g^\k) are coupling constants and are all possible operators constructed with 
the field <fi and n derivatives, which are compatible with the symmetries of the theory. The 
index i is used here to label different operators with the same number of derivatives. 

From the definition given above, it is easy to show that the functional satisfies the 
following "Exact Renormalization Group Equation" (or ERGE) Q] 



,dT k 1 
k— - = -Tr 

dk 2 



where the trace in the r.h.s. is a sum over the eigenvalues of the operator O (in flat space 
it would correspond to a momentum integration) and we have introduced the notation 
rj. 2 '* = tttt for the inverse propagator of the field defined by the functional IV The r.h.s. 



k 

of (Q can be regarded as the "beta functional" of the theory, giving the /c-dependence of 
all the couplings of the theory. In fact, taking the derivative of ([8]) one gets 

^EEa'"^"'- do) 

n=0 i 

where 

are the beta functions of the (generally dimensionful) couplings. Here we have introduced 
t = \og(k/ko), ko being an arbitrary initial value. If we expand the trace on the r.h.s. of 
in operators C\ and compare with (fit)]) , we can in principle read off the beta functions of 



the individual couplings. 

This formalism can be easily generalized to the case of multicomponent fields $ , where 
A may denote both internal and spacetime (Lorentz) indices. In this case the cutoff term will 
have the form <& A RkAB^ B , arid the trace in the r.h.s. of the ERGE will also involve a finite 
trace over the indices A, B. In the case of gauge theories there are further complications 
due to the fact that the cutoff interferes with gauge invariance. One way of dealing with 
this issue, which we shall use here, is to use the background field method [2^] (for another 
approach see 21[). One defines a functional of two fields, the background field and the 
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classical field, which is invariant under background gauge transformations. In the end the 
two fields are identified and one obtains a gauge invariant functional of one gauge field only. 

The ERGE can be seen formally as a RG-improved one loop equation. To see this, recall 
that given a bare action S (for a bosonic field), the one loop effective action is 



r« = S + ^Trlog 
Let us add to S the cutoff term ([7]); the functional 



5 2 S 



(12) 



r 



(i) 



S + -Trlo£ 



5 2 S 



+ Rk 



(13) 



may be called the "one loop effective average action" . It satisfies the equation 



(i) 



dk 



-Tr 



5 2 S 



Rk 



k 



dRk 
dk 



(14) 



which is formally identical to (Q except that in the r.h.s. the renormalized running cou- 
plings gi{k) are replaced everywhere by the "bare" couplings gi, appearing in S. Thus the 
"RG improvement" in the ERGE consists in replacing the bare couplings by the running 
renormalized couplings. In this connection, note that in general the cutoff function Rk may 
contain the couplings g^ and therefore the term k-^Rj. in the r.h.s. of will itself contain 
the beta functions. Thus, extracting the beta functions from the ERGE generally implies 
solving an algebraic equation where the beta functions appear on both sides 5 . This com- 
plication can be avoided by choosing the cutoff in such a way that it does not contain any 
coupling. Then, the entire content of the ERGE is in the (RG-improved) one loop beta 
functions. The result is still "exact" insofar as one is able to keep track of all possible 
couplings of the theory. 

The formal derivation of the ERGE from a path integral makes use of a bare action S. 
To explore the relation between S and the (renormalized) average effective action T^, for 
any k, requires that an ultraviolet regularization be defined (in addition to the infrared 
regularization provided by Rk). We will not need to discuss this point, since the bare action 
is an unphysical quantity and all the physics is encoded in the running renormalized action 
Tk- In this connection note that the trace in the r.h.s. of (Q, which includes an integration 

5 For example, in a scalar theory with action J d A x \Z(d(j)) 2 + m 2 <j> 2 + X(j) 4 ] it is natural to choose the cutoff 
of the form Rk(z) = Zkk 2 r(z/k 2 ). Then, the r.h.s. of (JHJ) will contain d t Z. 
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over momenta, is perfectly ultraviolet convergent and does not require any UV regulator. 
This is because the function k-^Rk in the r.h.s. goes to zero for momenta greater than k and 
makes the integration convergent. So, we can regard the derivation given above as merely 
formal manipulations that motivate the form of the ERGE, but then the ERGE itself is 
perfectly well defined, without the need of introducing an UV regulator. If we assume that 
at a given scale k physics is described by a renormalized action T^, the ERGE gives us a 
way of studying the dependence of this functional on k, and the behavior of the theory at 
high energy can be studied by taking the limit of for k — ► oo (which need not coincide 
with the bare action S). 

In most cases it is impossible to follow the flow of infinitely many couplings and a common 
procedure is to consider a truncation of the theory, namely to retain only a finite subset of 
terms in the effective action r*.. For example one could consider the derivative expansion 
(JHJ) and retain all terms up to some given order n. Whatever the choice, one calculates the 
coefficients of the retained operators in the r.h.s. of (Ej) and in this way the corresponding 
beta functions are computed. In general the set of couplings that one chooses in this way 
will not be closed under RG evolution, so one is neglecting the potential effect of the ex- 
cluded couplings on the ones that are retained. Still, in this way one can obtain genuine 
nonperturbative information, and this procedure has been applied to a variety of physical 
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problems, sometimes with good quantitative results. For reviews, see 

If we truncate the effective action in this way, there is usually no small parameter to 
allow us to estimate the error we are making. One indirect way to estimate the quality 
of a truncation relies on an analysis of the cutoff scheme dependence. The effective action 
Tfc obviously depends on the choice of the cutoff function R^. This dependence is similar 
to the scheme dependence of the renormalized effective action in perturbative QFT; only 
physically observable quantities derived from Tf, must be independent of Rk- This provides 
an indirect check on the quality of the truncation. For example, the critical exponents 
should be universal quantities and therefore cutoff-independent. In concrete calculations, 
usually involving a truncation of the action, critical exponents do depend on the cutoff 
scheme, and the observed dependence can be taken as a quantitative measure of the quality 
of the approximation. Ultimately, there is no substitute for performing calculations with 
truncations that contain more terms. Note that a good truncation is not necessarily one for 
which the new terms are small, but one for which the effect of the new terms on the old ones 
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is small. In other words, in search of a nontrivial FP, we want the addition of new terms 
not to affect too much the FP value of the "old" couplings, nor the "old" critical exponents. 



III. MATTER FIELDS AND CUTOFF SCHEMES 

In this section we illustrate the method that is used to compute the trace in the r.h.s. of 
dH]) in a gravitational setting and to evaluate the beta functions of the gravitational couplings. 
Quite generally, we will consider the contribution of fields whose inverse propagator is a 
differential operator of the form A = —V 2 + E, where V is a covariant derivative, both with 
respect to the gravitational field and possibly also with respect to other gauge connections 
coupled to the internal degrees of freedom of the field, and E is a linear map acting on 
the quantum field. In general, E could contain mass terms or terms linear in curvature. 
For example, in the case of a nonminimally coupled scalar, E = where £ is a coupling. 
A priori, nothing will be assumed about the gravitational action and also the spacetime 
dimension d can be left arbitrary at this stage. 

In order to write the ERGE we have to define the cutoff. For the operator to be used in 
the definition of ([7]), several possible choices suggest themselves. Let us split E = Ei + E 2 , 
where Ej does not contain any couplings and E 2 consists only of terms containing the 
couplings. We call a cutoff of type /, if Rk is a function of the "bare Laplacian" —V 2 , of type 
II ii it is a function of —V 2 + Ei and of type III if it is a function of the full kinetic operator 
A = —V 2 + E. The substantial difference between the first two types and the third is that 
in the latter case, due to the running of the couplings, the spectrum changes along the flow. 



For this reason these cutoffs are said to be "spectrally adjusted" [25|. 6 

Let us now restrict ourselves to the case when E 2 = 0, i.e. the kinetic operator does not 
depend on the couplings; then there is only a choice between cutoffs of type I and II. The 
derivation of the beta functions is technically simpler with a type II cutoff. In this case 
we choose a real function Rk with the properties listed in section II and define a modified 
inverse propagator 

P fc (A) = A + i4(A) . (15) 



6 In ||7j) it was assumed for simplicity that the operator O appearing in the argument of the cutoff function 
is also the operator whose eigenfunctions are used as a basis in the evaluation of the functional trace. It 
is worth stressing that this need not be the case, as discussed in Appendix A. 
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If the operator E does not contain couplings, using (1A10I) the trace in the r.h.s. of the ERGE 
reduces simply to: 

d t R k (A) 1 ~ (d t R k 



Tr- 



P fc (A) (Air)*/* 



^EQf.(f)^( A ) (16) 



j=0 

where B 2 i{A) are the heat kernel coefficients of the operator A and the Q-functionals, defined 
in ( 1A14IIA15|) are the analogs of momentum integrals in this curved spacetime setting. We 
have written dtRk to denote the derivative with respect to the explicit dependence of R k 
on k; when the argument of R k does not contain couplings this coincides with the total 
derivative ^R k - 

With a type I cutoff we use the same profile function R k but now with —V 2 as its 
argument. This implies the replacement of the inverse propagator A by 

A + R k (-V 2 ) = P k (-V 2 ) + E . (17) 

Therefore the r.h.s. of the ERGE will now contain the trace Tr pJ^^J^ . Since E is linear 
in curvature, in the limit when the components of the curvature tensor are uniformly much 
smaller than k 2 , we can expand 



P k + E f- 1 ' P[ +1 



e=o k 

Each one of the terms on the r.h.s. can then be evaluated in a way analogous to (lAlOj) . so 
in this case we get a double series: 

^-v^+'e = (i^sEEOH (p) / . (18) 

In order to extract the beta functions of the gravitational couplings one has to collect terms 
with the same monomials in curvature. We will see an example of this shortly. 

Before discussing specific examples, however, it is interesting to consider the scheme- 
independent part of the trace. In general, on dimensional grounds, the functionals Q n (^Mr^j 
appearing in ffTol) and (JTSJ) will be equal to k 2 ( n ~ m + 1 ) times a number depending on the profile 
function. As discussed in Appendix A, the integrals with m = n + 1 are independent of the 
shape of R k - Thus, in even-dimensional spacetimes with a cutoff of type II, and using ( 1A19I) . 
the coefficient of the term in the sum (TT6]) with % = § is Q (^pjf) B d (A) = 2B d (A). On 
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the other hand with a type I cutoff, using flAT8|) . flAT9|) and flM]) the terms with 
add up to 

d/2 



— I 



2 J dx y/gtr 
2B d (-V 2 + E 



b d (-V 2 ) - E6,_ 2 (-V 2 



1)d/2 -E d / 2 6 (-V 2 ^ 



{d/2) 



Therefore, in addition to being independent of the shape of the cutoff function, these coef- 
ficients are also the same using type I or type II cutoffs. 

As an example we will now specialize to four-dimensional gravity coupled to ns scalar 
fields, n e> Dirac fields, um gauge (Maxwell) fields, all massless and minimally coupled: 



d 4 Xy/g 



(19) 



In addition must contain a generic action for gravity of the form ([H]), which we do not 
write. (For the terms with four derivatives we use the parametrization given in (1761) below.) 
We shall compute here the contribution of these matter fields to the gravitational beta 
functions. The contribution of the gravitational field to its beta functions will be calculated 
in the next section using the same methods; as we shall see, the details of the calculation 
are technically more involved, but conceptually there is no difference. 

The field equation of each type of field defines a second order differential operator A^ = 
_ V 2 + E( A ), with A = S, D, M, gh and 



E (5) = 



E(D > = | ; 



E^ = Ricci 



(20) 



Here "Ricci" stands for the Ricci tensor regarded as a linear operator acting on vectors: 
Ricci(w) At = R^v v . For the gauge fields we have chosen the Lorentz gauge, and A^ 9 ^ is the 
operator acting on the scalar ghost. (It can be shown that the results do not depend on the 
choice of gauge 



With a type II cutoff, for each type of field we define the modified inverse propagator 
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P k (A^) = A^ + R k (AM). Then, the ERGE reduces simply to 



dt 



d t R k (A^)\ n D 



Pi 



(AW) J 



— Tr 



9 t E fc (A( c ))\ 



(AW) J 



2 

1 1 



^ ^ P fc (AW) J " nMTr(9h) V P t (AW) J 



(n 5 - An D + 2n M ) Q 2 



1 /OR 

+-R {n s + 2n D - 4n M ) Qi ± 5 h - 
6 \ ft 



180 



(3n s + 18n D + 36n M ) C 2 - (n s + llra D + 62n M ) £ 



+5n 5 P 2 + 12 (n s + n D - 3n M ) V 2 R 



(21) 



where C 2 is the square of Weyl's tensor and E = R pi/pa R pupa — AR^R^ + R (in four 
dimensions x = 32^2 / d 4 x ^/gE is Euler's topological invariant). The terms with zero and 
one power of R depend on the profile function R k but using (1A19I) we see that the coefficients 
of the four-derivative terms, i.e. the beta functions of gf\ are scheme-independent. 



With type I cutoffs, the modified inverse propagators are A^) + Rk(—V 2 ) 
E( A ) and the ERGE then becomes: 



ft(-V 2 ) + 



dTk 
dt 



n s 



Tr 



{S) 



dtP k (-V< 
Pk(-V 2 ) 



n D 



Tr 



(D) 



d t R k {-V 2 ) 
7 2 ) + f 



P k 



+ ^Tr 



(M) 



d t Rk\ 



- n M Ti [gh) 



d t Rk(-^ 



(22) 



v.p fc (-v 2 ) + Ricci; (9n> \p k (-v 2 ) 

Expanding each trace as in (lAlOp . collecting terms with the same number of derivatives of 
the metric, and keeping terms up to four derivatives we get 



dTk 
dt 



1 



2 (4tt) 



d 4 Xyfg 
d t Rk 



P k 

d t R k 
Pk 



(n s - An D + 2n M ) Q 2 

-Q2 



ns - { s Qi hr 



d t Rk 
Pk 

d t R k 

r k 



n D 



d t R k 
r k 



n M 



R 



1 

180 



(3n s + l8n D + 36n M ) C 2 - (n s + Un D + 62n M ) E 
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+5n s R 2 + 12 (n s + n D - 3n M ) V 2 i? 



(23) 



We see that the terms linear in curvature, which contribute to the beta function of New- 
ton's constant, have changed. However, the terms quadratic in curvature have the same 
coefficients as before, confirming that the beta functions of the dimensionless couplings are 
scheme-independent . 

In order to have more explicit formulae, and in numerical work, one needs to calculate 
also the scheme-dependent Q-functionals. This requires fixing the profile Rk- In this paper 
we will mostly use the so-called optimized cutoff flA2ip in which the integrals are readily 
evaluated, see equations (!A22fA23|IA24p . This cutoff has the very convenient property that 
Q_ n (^) = for n > 1. Thus, the sum over heat kernel coefficients on the r.h.s. of (lAlOj) 



\ Pk J 

terminates. In particular, in four dimensions, there are no terms beyond those that are 
explicitly written in (pZTj) or (|23|) 7 . For more general cutoffs a calculation of beta functions 
for curvature-polynomials of cubic and higher order would require the knowledge of higher 
heat kernel coefficients. 

Let us briefly comment on the spectrally adjusted (type III) cutoffs. These only occur 
when the kinetic operator A contains couplings, for example a mass term or, in the case of 
a scalar field, a nonminimal coupling of the form £R. In this case the last factor in (Q is 

^ = d t R k + Y,R' k ^d tg ,, 

where R' k denotes the derivative of Rk{z) with respect to z and the sum extends over all 
couplings (this assumes that the derivative of the operator appearing in Rk commutes with 
the operator itself). This introduces further nonlinearities in the system. Since the beta 
functions of the couplings appear linearly in the r.h.s. of the equation, to obtain the beta 
functions one has to solve a system of linear equations. 

In a consistent truncation one would have to add to the terms in (1231 or fl2Tj) the con- 
tribution due to the gravitational field. This will be done in the next section. For the time 
being we observe that if the number of matter fields is of order A" — > oo, this is the dominant 
contribution and constitutes the leading order of a 1/N expansion 



The matter con- 



tributions by themselves have a form that leads to a gravitational FP. Comparing equation 



' this had also been observed in a different context in 



27] 



1(3 



(jTOj) with equations (123|) or (|2T|) one can read off the beta functions, which all have the form 

a^k 4 



A (») 

d 9i An),4- 



dt 

where a± are constants. Then, the beta functions of the dimensionless variables = 

k n ~ 4 gl n ^ are 

\ = {n- + «?° . (24) 

This simple flow has indeed a FP for all couplings. For n^4 

(n) 

= ^ , (25) 
in particular writing g(°' = 2ZA and g^ = —Z = —ttt-^, the flow of the cosmological 



constant and Newton's constant is given by 

d~A 



-2A + 87ra (0) G + 167ra (2) GA , 
dt 

^ =2G + 167ta^G 2 , (26) 
dt 

which, for a type II cutoff, has a FP at 

K = J_ns-±n D + 2n M ^ = 12tt 

4% + 2n D - An M —n s - 2n D + An M 

Note that the FP occurs for positive or negative A depending on whether there are more 
bosonic of fermionic degrees of freedom. The FP value of G, on the other hand, will be 
positive provided there are not too many scalar fields. For n = 4, (T2"4"l) gives a logarithmic 
running 

g i t\k)=gt\h) + af ) Hk/k Q ) , 

implying asymptotic freedom for the couplings l/gf^. This is the same behavior that is 
observed in Yang-Mills theories and is in accordance with earlier perturbative calculations 
0, [3^]. As noted, it follows from (1A24|) that with the optimized cutoff, for n > 4, g^ = 0. 



The critical exponents at the nontrivial FP are equal to the canonical dimensions of the 
g^'s, so A and G are UV-relevant (attractive), l/g^ are marginal and all the higher terms 
are UV-irrelevant. Note that in perturbation theory G is irrelevant. At the nontrivial FP 
the quantum corrections conspire with the classical dimensions of A and G to reconstruct 
the dimensions of g(°> and g( 2 >. This must happen because the critical exponents for g(°> and 
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FIG. 1: The generic form of the flow induced by matter fields. 



g^> are equal to their canonical dimensions and the critical exponents are invariant under 
regular coordinate transformations in the space of all couplings; the trasformation between 
G and is regular at the nontrivial FP, but it is singular at the Gaufiian FP. 

This simple flow is exact in the limit iV — * oo, but is also a rough approximation when 
graviton effects are taken into account, as we shall discuss in sections IV-G and VI. It is 
shown in figure [TJ 



IV. EINSTEIN'S THEORY 



As a first step towards the inclusion of quantum gravitational effects, we discuss in this 
section the RG flow for Einstein's gravity, with or without cosmological constant. This 
truncation has been extensively discussed before 



10j. Here we will extend those results 



in various directions. Since the dependence o 
function Rk has already been discussed in 



the results on the choice of gauge and profile 



10 



111 ] here we shall fix our attention on a 



particular gauge and profile function, and analyze instead the dependence of the results on 
different ways of implementing the cutoff procedure. The simplicity of the truncation will 
allow us to compare the results of different approximations and cutoff schemes, a luxury 
that is progressively reduced going to more complicated truncations. 

The theory is parametrized by the cosmological constant A and Newton's constant G = 
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1/ (16nZ), so that we set = 2AZ and gW = —Z in equation (jSJ). All higher couplings 
are neglected. Then the truncation takes the form 

r fc = J dx^ (2AZ - ZR(g)) + S GF + S ghost , (28) 

where Sgf is a gauge-fixing term and Sghost is the ghost action. We decompose the metric 
into g pv = gl„ ■+h /ll/ where Qiiv is a background. We will refer to the field as the graviton, 
even though it is not assumed to be a small perturbation. We consider background gauges 
of the type: 

S GF (g {B) ,h) = A j dx^x^g^x, , (29) 



where 

1 + p 
d 



All covariant derivatives are with respect to the background metric. In the following all 
metrics will be background metrics, and we will omit the superscript (B) for notational 
simplicity. In this section we will restrict ourselves to the de Donder gauge with parameters 
a = 1, p = | — 1, which leads to considerable simplification. The inverse propagator of h^, 
including the gauge fixing term, can be written in the form 

i J dx ^/gh pv T { ^ upa h pa 

containing the minimal operator: 

rl 2) ;; = ^[^(-v 2 -2A) + f/;;] , (30) 

where 8 

U% = RK% + \ (g^R pa + R^g pa ) - #J£> - R^ a) . 

In the following we will sometimes suppress indices for notational clarity; we will use boldface 
symbols to indicate linear operators on the space of symmetric tensors. For example, the 
objects defined above will be denoted K, 1, P, U. Note that P and 1 — P are projectors 
onto the trace and tracefree parts in the space of symmetric tensors: h pv = h^f^ + hffi 



These definitions coincide with those of [8| except that the A term has been removed from U. 
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where h$) = P^h pa = jg^h. Using that K = \ ((1 - P) + ^P), if d ^ 2 we can rewrite 
equation (1301 in either of the following forms: 



f ] = ZK(-V 2 -2A1 + W) 
Z 



1 - P) (-V 2 - 2A1 + 2U) 



d-2. 



-V 2 - 2A1 - 



d-2 



U 



(31) 



where we have defined 



W% = 2U% - ^L^L{R pa g^ + g^RT - Rg pa9 n . 

Note that the overall sign of the second term in the second line of (13T]) is negative when 
d > 2. This is the famous problem of the unboundedness of the Euclidean Einstein-Hilbert 
action. We shall see shortly how this is dealt with in the ERGE. Later on, we will need the 
traces: 



tr 1 



d(d+l] 



trP = 1 : tr (1 - PI 



d 2 + d - 2 



; trW 



d(d- 1) 



R; 



trW 2 = 3iWi^ + * ~ M + + * ~ ¥ + M + 4 & 



d-2 



2{d-2) 



(32) 



The ghost action is: 



^ g host = - J s/g~c, (-v 2 ^ - K) c v . 

On the c?-dimensional sphere we can write 

1 



U 



,d 2 -3d + 4 rf-2d-4„ 

1 - P) ~ — ; —R ~ P —R 



d(d-l) 



d 



Then, using the second line of (13 ip . we have 



,(2) 



1 - P) -V 2 - 2A 



d 2 - 3d + 4 
d(d-l) 



R 



d-2. 



2A 



d-4 



d 



-R 



(33) 



We will now discuss separately various types of cutoff schemes. 



A. Cutoff of type la 



This is the scheme that was used originally in |8j. It is defined by the cutoff term 
AS^Vl = \J dx^h pu R k (-V 2 r upa V ~ Jdx s/gC^i-^YvC" , (34) 
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where 



R 



(9h), 



-v 2 r, 



ZKR k (-V 2 ) 



(35) 



for gravitons and ghosts respectively. Defining the anomalous dimension by 

1 dZ 



71 Z dt 



we then have 



dUk 
~df 



ZK [d t R k (-V 2 )+r)R k (-V< 



(36) 



(37) 



The calculation in |8j proceeded as follows. The background metric is chosen to be that 
of Euclidean de Sitter space. The modified inverse propagator is obtained from (|33|) just 
replacing —V 2 by P k (— V 2 ). Using the properties of the projectors, its inversion is trivial: 



-.(2) 



R; 



2 
Z 



P k -2k+^^R d-2 Pk -2A+^R 



(38) 



Decomposing in the same way the term j|Rfc, multiplying and tracing over spacetime indices 



one obtains 

dT k _ 1 
dt 2 



Trfl 



d t R k + r]R k 



l TrP d t R k + V R k _ m_ 
2 P k -2A + ^R u P k 



k d - K d 

One can now expand to first order in R, use the traces (132]) and formula flAlOj) to obtain: 



dT k 
dt 



( 4vr )d/2 

d (d 



dx^/g< 



d(d 



+ 



24 
d(d-l] 



Q d _ 



4 

dtRk + yRh 
Pk - 2A 
d t R k + r]R k 



1)^ fd t R k + V R k 
~ Q n R-2A 



— Q d i 

6^2"! 



-dQ 
dtRk 



d t Rk 
P k 



Q 



dtRk 



R + 0{R 2 



(39) 



(Pk ~ 2A)' 

This derivation highlights two noteworthy facts. The first is that the negative sign of the 
kinetic term for the trace part of h is immaterial. With the chosen form for the cutoff, 
any prefactor multiplying the kinetic operator in the inverse propagator cancels out between 
the two factors in the r.h.s. of the ERGE. The second fact, which we will exploit in the 
following, is that the singularity occurring in the kinetic operator for the trace part in d = 2 
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(see equation (!3T|) ) is actually made harmless by a hidden factor d — 2 occurring in U. So, 
the final result fl39l) is perfectly well defined also in two dimensions. 

On the other hand, an issue that is sometimes raised in connection with this calculation 
is background dependence. The calculations in section III were done without choosing a 
specific background, so the question arises whether the same can be done here. The answer 
is positive, provided we do not decompose the field h^ v into tracefree and trace parts, and 
we use for the inverse propagator the form given in the first line of (l3~Tj) . Then, the modified 
inverse propagator for gravitons is 



Tf + R k = ZK (P fc (-V 2 ) - 2A1 + W) 



(40) 



(2) 

On a general background it is impossible to invert T k + exactly, but remembering that 
W is linear in curvature we can expand to first order: 



,(2) 



R; 



F (2) /i 
1 CC v 



-1 



-1 



Pk 
1 

Pk 



2A 



P k -2A 



W + 0(R 2 



-Tfc 



(41) 



Then the ERGE becomes, up to terms of higher order in curvature, 



dT k I d t R k + rjRk 
dt 2 1 P k - 2A 



1 



W 



Tr 



dtRk 



Pk 



St 1 , 



— i? M 

Pk 1 



P k -1t 

From here, using (IA10[) one arrives again at This alternative derivation explicitly 

highlights the background independence of the results. 

We are now ready to extract the beta functions. The first line of (1391) gives the beta 
function of 2ZA, while the other two lines give the beta function of —Z. Note the appearance 
of the beta function of Z in the r] terms on the r.h.s. In a perturbative one-loop calculation 
such terms would be absent; they are a result of the "renormalization group improvement" 
implicit in the ERGE. The beta functions can be written in the form 

k d 



d_ / 2A \ 
di [mrGJ 
d ( 1 



167T 
k d ~ 2 



{A 1 + A 2 v) 
(B 1 + B 2V ) 



(42) 



dt \16nG J 16vr 

where Ai, A 2 , R>i and B 2 are dimensionless functions of A, k and of d which, by dimensional 
analysis, can also be written as functions of A = Afc -2 and d. One can solve these equations 
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for ^ and obtaining 



dA _ _ 2K + ^A 1 + 2B 1 A + G(A 1 B 2 -A 2 B 1 ] 



dt 2(1 + B 2 G) 

d 4 = (d-2)G + ^L. (43) 
dt 1 + B 2 G 

The corresponding perturbative one loop beta functions are obtained by neglecting the r\ 
terms in (142]) . i.e. setting A 2 = B 2 = 0, and expanding A\ and B\ in A. The leading term is 

^ = -2A + ^4 1 (0)G + 5 1 (0)GA, 
dC 

_ = (d - 2)G + B^G 2 , (44) 

where A± and B\ are evaluated at A = 0. This flow has the same structure of the one written 
in (126]) ; we will refer to it as the "perturbative Einstein-Hilbert flow". We will discuss its 
solution in section IV G. 

The explicit form of the coefficients appearing in (1421) . with the optimized cutoff, is 

16tt(c/-3 + 8A) 

(47r)*r(f)(l-2A) 
16vr (d+1) 



A 



B l 



B 



(47r)'(d+2)r(f)(l-2A) 

-4tt(-c/ 3 + 15c/ 2 - 12d + 48 + (2d 3 - Ud 2 - 192)A + (16c/ 2 + 192)A 2 ) 

3(47r)^ (1 - 2 A) 2 

4tt (cl 2 - 9 d + 14 - 2 (rf + 1) (d + 2) A) 
3(47r)^(d + 2)r(f)(l-2A) 2 



A similar form of the beta functions had been given in 3JJ in another gauge. For the sake 
of clarity we write here the beta functions in four dimensions: 

G 3 - 4A - 12A 2 - 56A 3 + 

LZ7T 



Pa = -2A + 



6tt (i _ 2A) 2 - 

_ ^ & 11-18A + 28A 2 

&r = 2G = =— ■ . (45) 

3tt (i _ 2A) 2 - 

Note the nontrivial denominators, which in a series expansion could be seen as resummations 
of infinitely many terms of perturbation theory. They are the result of the "RG improve- 
ment" in the ERGE. 
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B. Cutoff" of type lb 



This type of cutoff was introduced in 32J. The fluctuation and the ghosts are de- 
composed into their different spin components according to 



h 



1 ^2 1 , 



and 



+ V"c ,C U 



(46) 



(47) 



where h/E v is tranverse and traceless, £ is a transverse vector, a and h are scalars, c T and 
c T are transverse vectors, and c and c are scalars. These fields are subject to the following 
differential constraints: 



h T » = 0; V v h 



; 



V v i v = ; Wc T u = ; V„c T ^ = . 



Using this decomposition can be advantageous in some cases because it can lead to a partial 
diagonalization of the kinetic operator and it allows an exact inversion. This is the case 
for example when the background is a maximally symmetric metric. In this section we 
will therefore assume that the background is a sphere; this is enough to extract exactly 
and unambiguously the beta functions of the cosmological constant and Newton's constant. 
Then the ERGE can be written down for arbitary gauge a and p. We refer to 32j for 
more details of the calculation. In the gauge a = 1 and without making any approximation, 
the inverse propagators of the individual components are 



r (2) 



p(2) 



Z 
~2 



-V 2 



d 2 - 3d + 4 



d(d 



Z -V 2 



R 
1 



■1) 
-V 2 



R-2A 
d-3 



d 



R-2A 



,(2) 
hh 



p(2) 



,(2) 



,(2) 



-z 



-id 



-V 2 + ^ 



.d 



2d 
7 + ~d 



(-v 2 ) 



d 

-V 2 - 



2A 



R 



d-4 



R-2A 



fj 



(48) 



The change of variables (146]) and (1471) leads to Jacobian determinants involving the operators 

R 



o R 



d 



1 J c 
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for the vector, scalar and ghost parts. The inverse propagators (|48j) contain four derivative 



terms. In 



10 



32] this was avoided by making the field redefinitions 



-V 2 



R 



-a. 



(49) 



d- 1 

At the same time, such redefinitions also eliminate the Jacobians. These field redefinitions 
work well for truncations containing up to two powers of curvature, but cause poles for 
higher truncations as the heat kernel expansion will involve derivatives of the trace argu- 
ments. Therefore, in later sections we will not perform the field redefinitions, but treat the 
Jacobians instead as further contribution to the ERGE by exponentiating them, introducing 
appropriate auxiliary fields and a cutoff on these variables. Here we describe the result of 
performing the field redefinitions. The ERGE is 

dT k 1^ d t R k + 7]R k 1 d t R k + r]R k 



dt 



2A 



d 2 -3d+4 p 
d(d-l) - rt 

d t R k + rjR k 



X Tr' 
2 ir « 



+ 2 Tr(0) P,-2A+^i? 



ZTr" 



d t R k 



2\+^R 



7]R k 



-Tr (1) 



k 

dtR k 



(o) 



2A + ^R 



Pk- 



R 



dtRk 



(0) p _ 2R 



(50) 



k d rk d 

The first term comes from the spin-2, transverse traceless components, the second from 
the spin-1 transverse vector, the third and fourth from the scalars h and a. The last two 
contributions come from the transverse and longitudinal components of the ghosts. A prime 
or a double prime indicate that the first or the first and second eigenvalues have to be 
omitted from the trace. The reason for this is explained in Appendix B. 

Expanding the denominators to first order in R, but keeping the exact dependence on A 
as in the case of a type la cutoff, and using the formula flAlOl) . one obtains 

dtRk 



d£k 
dt 



d A -2d 3 -d 2 -4d + 2 



+R 
d l 



+ 



Ad(d- 
13d 2 - 24d ■ 



1) 
12 



Q 



yRk 
Pk -2A 

d t R k + r/Rk 



dQd 

2 
d 



d t Rk 
Pk 



Qi-x 



(Pk - 2A) 2 
d t R k + vRk\ d 2 



j — Q± 

d 2 



Pk — 2A 



6 



6d 



Qi-i 



d t Rk 

Pi 

d t Rk 
Pk 



(51) 



+ 0(R 2 



24d(d- 1) 

In principle in two dimensions one has to subtract the contributions of some excluded modes. 
However, using the results in Appendix B, the contributions of these isolated modes turn 
out to cancel. Thus, the ERGE is continuous in the dimension also at d — 2. 



25 



The beta functions have again the form (1431) ; the coefficients A\ and A% are the same as 
for the type la cutoff but now the coefficients B\ and B 2 are 

B v = An (d(d - l)(d 3 - 15d 2 - 36) + 24 -2(d 5 -8d 4 -5d 3 -72d 2 -36d + 96) A 
-16(d - l)(d 3 + 6d + 12)A 2 ) /3(4n)ld 2 (d - 1)F (-) (f - 2A) 2 



Bo = 4tt 



2 y 

d{d 4 - 10d 3 + lld 2 - 38d + f 2) - 2{d + 2){d A -13 d 2 - 24d + 12) A 
3(4tt)! (2 + d)(d - l)d 2 r(f )(1 - 2A) 2 



In four dimensions, the beta functions are 

1 (12 - 33A + 20A 2 - 200A 3 )G + ^~^ G 2 

g. = _2A H - ' _ 12n 

24vr (1 _ 2A) 2 - ^G 

= 26 _^(105-212A + 200A 2 )G 2 " 
G 24vr (1 _ 2A) 2 - M<3 

In order to appreciate the numerical differences between this procedure and the one where 
the fields £ M and a are not redefined as in (1491) . we report in Appendix D the results of the 
alternative formulation. 



C. Cutoff of type II 

Let us define the following operators acting on gravitons and on ghosts: 

A 2 = -V 2 + W (53) 
A ( 9 h) = -V 2 " Ricci • (54) 
The traces of the b 2 -coefficients of the heat-kernel expansion for these operators are 

tab , (A> , = *(fi-w) = !SLiM* 



trb 2 (A g/l ) = trf-l + Riccij = —jr~R 



The type II cutoff is defined by the choice 



R fc = ZKR k (A 2 ) 
R k 9h M v = 5^R k (A( gh )) 



26 



which results in 



rf+Rfc = ZK(P,(A 2 )-2A) 
T% + Ri 9h) = P k (A (gh) ) (55) 



ZK(d t R k (A 2 )+rjR k (A 2 )) . 



and 

dH k 
~df 

Collecting all terms and evaluating the traces leads to 
dT k _ l Tr <9A(A 2 ) + yR k (A 2 ) ^d t R k (A igh) ) 



dt 2 P fc (A 2 )-2A Pk(A {gh) ) 



1 f jdid+T) ft 



d(7-5d)_ fd t R k + r!R k \ d + 6„ ( d t R k 

^#-1 D OA ^—^1-1 



P k -2A J ^ V ^ 

24 -f-V P.-2A y 6 .^ VPfc;j ^ + ^)}.(56) 

The beta functions are again of the form (I43p . and the coefficients Ai and 7^2 are the same 
as in the case of the cutoffs of type I. The coefficients B\ and B 2 are now 

4vr(5d 2 - 3d + 24 - 8(d + 6)A) 
1 " 3(4vr)lr(f)(l-2A) 
B = A7r(5d - 7) 

3(4vr)fr(f)(f-2A) 

In four dimensions, the beta functions are 

I (3 - 28A + 84A 2 - 80A 3 )G + 191 ~ 5 7r 12A & 



Px = -2A-r 
A 6vr (i_2A)(l-2A-^G) 

I (23 - 20A)G 2 

3^(l-2A)-^G 



D. Cutoff of type III 



(2) 

cutoff function as a function of the whole inverse propagator T k , only stripped of the overall 
wave function renormalization constants. In the case of the graviton, rj; 2) = ZK(A 2 - 2A1) 
while for the ghosts T^L, = A gh , where A 2 and A gh were defined in (1531) . Type III cutoff is 
defined by the choice 

R fc = ZKR k (A 2 - 2A) (58) 
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for gravitons, while for ghosts it is the same as in the case of type II cutoff. Since the 
operator in the graviton cutoff now contains the coupling A, the derivative of the graviton 
cutoff now involves an additional term: 



dH k 



ZK {d t R k (A 2 - 2A) + r]R k (A 2 - 2A) - 2R' k (A 2 - 2A)d t A) (59) 



dt 

where R' k denotes the partial derivative of Rk{z) with respect to z. Note that the use of the 
chain rule in the last term is only legitimate if the t-derivative of the operator appearing 
as the argument of R k commutes with the operator itself. This is the case for the operator 
A 2 — 2A, since its t-derivative is proportional to the identity. The modified inverse propagator 
is then simply 

r| 2) +R fc = ZKP fc (A 2 -2A) 
for gravitons, while for ghosts it is again given by equation (1551) . Collecting, 



dY k _ 1 d t R k (A 2 - 2A) + v R k (A 2 - 2A) - 2R' k (A 2 - 2A)d t A d t R k {A {gh) ) 
dt 2 1 P k (A 2 - 2A) 1 P k (A (gh) ) • 1 j 

The traces over the ghosts are exactly as in the case of a cutoff of type II. As in previous 

cases, one should now proceed to evaluate the trace over the tensors using equation (lAlOh 



and the heat kernel coefficients of the operator A 2 — 2A. However, the situation is now more 
complicated because the heat kernel coefficients B 2k {A 2 — 2 A) contain terms proportional 
to A fc and A h ~ l R, all of which contribute to the beta functions of 2AZ and —Z. This 
is in contrast to the calculations with cutoffs of type I and II, where only the first two 
heat kernel coefficients contributed to the beta functions of 2AZ and —Z. In order to 
resum all these contributions, one can proceed as follows. We define the function W(z) = 
a t R k ^)+vRk(z)-2R' k (z)d t A and the function w{ z ) = W(z-2A). It is shown explicitly in the end 

of Appendix A (equation ( 1A33I) and following) that TrW = TrW. Then, the terms without 
R and the terms linear in R (which give the beta functions of 2AZ and —Z respectively) 
correspond to the first two lines in (1A34I) . In this way we obtain 

dT k If, n /d(d+l)^(2AV fd t R k + V R k -2d t AR' k \ ^ fd t R 



1 f, r J (2A)* f d t R k + V R k - 2d t AR' k \ 



k 



d(7-M) ^(2A)' ( d t R t + r/Rt - 2d,AR' k \ rf + 6 (a,R„\ 1 

The remarkable property of the optimized cutoff is that in even dimensions the sums in those 
expressions contain only a finite number of terms; in odd dimensions the sum is infinite but 
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can still be evaluated analytically. Using the results (IA22IA23IA24IA29IA30IA31IA32!) the 
first sum in (l6Tj) gives 

1 d + 1 (k 2 + 2A) d / 2 / r] fc 2 + 2A , n d t A 
^ir) d / 2 ~2 V{dj2) I + |TT W~ 

whereas the second sum gives 



2 J^)J dx ^ (62) 



1 d(7-5d) (fc 2 + 2A)^ / r7 fc 2 + 2A d t A\ 

/2 + ^^^ + 2^ e^V^ (63) 



(4 7r )d/2 24 r(d/2) V d / 2 fc2 k2 

This resummation can actually be done also with other cutoff's. An alternative derivation 
of these formulae, based on the proper time form of the ERGE is given in Appendix C. 

The beta functions cannot be written in the form (I43p anymore, because of the presence 
of the derivatives of A on the right hand side of the ERGE. Instead of (142 j) we have 

(A 1 + A 2 rj + A 3 d t A) , 



dt \16nG J 16tt 
d f 1 \ k 
"dt \ IQtiG 



(lib) + (64) 



where 



A, 



167r(-4+(rf+l)(l + 2A)2 +1 ) 

(4vr)fr(f) 
167r(d+l)(l + 2A)f +1 

(47r)^(d + 2)r(f) 
167r(d+l)(l + 2A)3 

47r(-4(rf+ 6) + d(7 - 5rf)(l + 2A)i; 
3(4vr)fr(f) 

47T(7-5rf)(l + 2A)^ 

3(47r)#r(f) 
47rrf(7-5rf)(l + 2A)i- 1 



3(4vr)lr(f) 

Solving flMl) for dX/dt and dG/dt gives 

rfA = _ 2A (Ai + 2(gi - AQA - 4B 3 X 2 )G + (M - Aggi + 2(A 2 B 3 - A 3 B 2 )X)G 2 

dt ' 2 + (2B 2 - A 3 - 2B 3 X)G + (A 2 B 3 - A 3 B 2 )G 2 

dG _ & | 2(B 1 - 2B 3 X)G 2 + {MB, - A 3 B 1 )G 3 
dt 2 + (2B 2 -A 3 - 2B 3 X)G + {A 2 B 3 - A 3 B 2 )G 2 ' 
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Cutoff type lb, d=4 




FIG. 2: The beta function of G with A = and cutoffs of type la and lb. The perturbative one 
loop result in light gray, the RG-improved one in darker color. 



In four dimensions, the beta functions are 



1 (3 + 14A + 8A 2 )G + fi±|£t(i9i - 60A - 260A 2 )G 
6^ l- I ^(43 + 120A + 68A 2 )G + 7 f^(l + 2A) 4 G 2 



-2A + 



(23 + 26A)G 2 - 



2 51+152A+100A' 



G 3 



3vr l _ ^(43 + 120A + 68A 2 )G + ^(1 + 2A) 4 G 2 



(66) 



E. Gravity without cosmological constant 



Before discussing the general case, it is instructive to consider the case A = 0. In terms 
of the dimensionless coupling G the beta functions with cutoffs of type I and II have the 
form (l4"3"j) . where the constants B\ and B 2 are evaluated at A = 0. The cutoff of type III 
leads instead to the more complicated beta function (1651) . with A = 0. The beta function 
of G in four dimensions is shown in figures [2] and [3] (blue, dark lines) for different cutoff 
types. It always has a Gaufiian FP in the origin and a nontrivial FP at G* = — Bi+ ^ 2 )b 2 ' 
The Gaufiian FP is always UV-repulsive (positive slope) whereas the non-Gaufiian FP is 
UV-attractive. In a theory with a single coupling constant such as this, the slope of the 
beta function at the nontrivial FP is related to u, the mass critical exponent. For type I 
and II cutoffs it is given by 



v 



d(3, 



G 



dG 



(d-2)(l + (d-2)^ 



(67) 



The leading, classical term is universally equal to d — 2, the correction is scheme-dependent. 
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Cutoff type II, d=4 C Cutoff type III, d=4 




FIG. 3: The beta function of G with A = and cutoffs of type II and III. The perturbative one 
loop result in light gray, the RG-improved one in darker color. In the case of cutoff type III the 
vertical line is the asymptote of the RG-improved beta function. 

It is instructive to compare these "RG improved" results to the respective "one loop 
approximations". As discussed in section 2, the "unimproved" perturbative one-loop beta 
functions are obtained by neglecting the derivatives of the couplings occurring in the r.h.s. 
of the ERGE, in practice setting B 2 = 0: 

08\uoop= {d-2)G+B 1 G 2 . (68) 

These beta functions are shown as the gray, light lines in figures |2] and [3J Since B\ < 
0, they are inverted parabolas, with a GauBian FP in the origin and a nontrivial FP at 

= — (d — 2)/Bi. The slope at the nontrivial FP is always the opposite of the one at the 
GauBian FP, and therefore equal to —2. Note that both (j3j) and ([6]) are of this form, for 
specific values of the constant B\. Because B\ < B 2 < 0, the "RG improved" FP occurs 
always at smaller values of G* than the corresponding perturbative one. The GauBian FP 
is always UV-repulsive (positive slope) whereas the non-GauBian FP is UV-attractive. In 
table I we report the numerical values of the coefficients Bi and B 2 and the position of the 
FP and the critical exponent in four dimensions. 

One can see from figure [2] that for cutoffs of type I the one-loop approximation is quite 
good up to the nontrivial FP and a little beyond. For larger values of G the effect of the 
denominator in ( 143|) becomes important; the beta function deviates strongly from the one- 
loop approximation and has a negative pole at G = — 1/B 2 > G*. When one considers type 
II and III cutoffs, the RG improved beta functions deviate from the perturbative one loop 
beta functions sooner and the effect is stronger; in the case of the type III cutoff the pole 
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Approximation 


Bi(0) 


B 2 (0) 


G* 




e exp., leading order 






0.158 


2.000 


la 

la - 1 loop 


n 

OvT 


l 

1 9-TT 


1.639 
1.714 


2.091 
2.000 


lb with f.r. 

Ib with f.r. - 1 loop 


35 

O/T 


29 
72tv 


1.213 
1.436 


2.368 
2.000 


lb without f.r. 

Ib without f.r. - 1 loop 


131 

457T 


61 

79n-7r 

i ZU7T 


2.040 
2.158 


2.116 
2.000 


II 

II - 1 loop 


23 

37T 


13 

127T 


0.954 
0.820 


2.565 
2.000 


III 

III - 1 loop 


23 
3tt 


13 
12?r 


0.424 
0.820 


3.870 
2.000 



TABLE I: The fixed point for Einstein's theory in d = 4 without cosmological constant. The 
leading beta function for the e expansion is derived in section IVF; the one for the cutoff of type 
Ib without field redefinitions is given in Appendix D. 

occurs before the one loop beta function has the zero. We see that the RG improvement 
leads to stronger effects if more terms of the operator are taken into the definition of the 
cutoff. It is interesting to observe that the RG improvement always brings the nontrivial 
fixed point closer to the perturbative regime. Since at low energies G is close to zero, the 
region of physical interest is < G < G*. Thus, the pole in the beta functions at finite 
values of G should not worry us. 

F. The e expansion. 

Before discussing the four dimensional case, it is useful and instructive to consider the 
situation in arbitrary dimensions. In particular, this will allow us to compare the results of 
the ERGE with those obtained in the e-expansion. We have seen within the approximations 
of section 2 that in d dimensions the beta function of the dimensionless coefficient of the R d ^ 2 
term is scheme-independent. Therefore in two dimensions one expects the beta function of 
Newton's constant, or at least its leading term, to be scheme-independent. This is confirmed 
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FIG. 4: The position of G> as a function of d in the one loop approximation with cutoff of type la 
(central line in black). The other lines are the first five orders of the e expansion (order n means 
the beta function has been expanded to order n in e.) The zeroth order is the light gray straight 
line. Higher orders are represented by darker shades of gray. Order 1 and order 3 have singularities 
at d » 2.998 and d « 3.609. 



ft|«=~ (») 



by formula (168]) for the one loop beta function, and the results listed in sections 4.1 to 4.4: 

38 
T 

for all types of cutoff. On the other hand the coefficient B2 is scheme-dependent. This 
mirrors the well known fact that in perturbation theory the leading term of the beta functions 
of dimensionless couplings is scheme-independent and higher loop corrections are not. 

The beta function (143!) with A = can be solved exactly for any d. The nontrivial FP 
occurs at: 

* = - Bt + V-QB, ■ (70 » 
Knowing the solution in any dimension we can now check a posteriori how good the e- 

expansion is. To this end we have to expand the beta function in powers of e = d — 2 and 
look for FPs of the approximated beta functions. 

The leading term of the e expansion consists in retaining only the constant, scheme- 
independent term fl69|) . Then the beta function is given by equation (jSJ) and the fixed point 
occurs (for any cutoff type) at 
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Cutoff type la 



G, 




Cutoff type lb 



0.05 0.1 



0.15 0.2 



0.25 




FIG. 5: The position of the nontrivial fixed point as a function of d for cutoffs of type la (left 
panel) and lb (right panel). 

We see from table I that this value is quite small compared to the direct calculations in 
four dimensions. For the higher orders of this expansion it is necessary to specialize the 
discussion to a specific type of cutoff. For definiteness we will consider the case of a type la 
cutoff, for which (setting A to zero) 

4tt (d 3 - 15c/ 2 + 12d - 48) 



To order e the beta function is 
38 ~, 



G 



"T G2 + 



G 



3(47r) d / 2 cir(f) 



1 197 + 381og2 t 191ogvr 



(71) 



G 



e + 0(e 5 



3 3 3 3 

Whereas to leading order the solution exists for all d, to this order the solution has a 
singularity for e « 2.99679. The occurrence of such singularities at finite values of e is 
expected. The one loop solution together with some of its approximants is plotted in figure 
H] for dimensions 2 < d < 4. We see that when the beta function is expanded to even order 
in e (to = 0, 2, 4 in the figure) the e expansion significantly underestimates the value of G* 
in d — 4 whereas for odd order (to = 1, 3 in the figure) it has a positive pole at some value 
of the dimension. On the other hand from equations ( 17D1) and ( 17TT) one sees that for large d, 
G* grows faster than exponentially. 

We conclude this discussion by mentioning that when the e expansion is used in presence 
of a cosmological constant, there are several FPs and even for e very small they have negative 
A*. Thus the e expansion is not very helpful in the presence of A. One can solve exactly the 
equations 4A = and ^tG = for arbitrary d and plot the position of the FP in the A-G 
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Cutoff type II G - Cutoff type III 




FIG. 6: The position of the nontrivial fixed point as a function of d for cutoffs of type II (left 
panel) and III (right panel). 

plane as a function of d. This is shown in figures [5] and [61 The fixed point is in the origin 
at d = 2; as d grows, G* grows monotonically while A* is initially negative, then becomes 
positive. For moderately large dimensions (of order 10) becomes very large (of the order 
10 6 ) while A* < 1/2 always. 



G. Four dimensions 



Let us now consider Einstein's theory with cosmological constant in d = 4. The beta 
functions for A and G for various cutoff types have been given in equations (I45ll52|57ll66l) . 
All of these beta functions admit a trivial (Gaufiian) FP at A = and G = and a 
nontrivial FP at positive values of A and G. Let us discuss the Gaufiian FP first. As usual, 
the perturbative critical exponents are equal to 2 and —2, the canonical mass dimensions 
of A and G. However, the corresponding eigenvectors are not aligned with the A and G 
axes. It is instructive to trace the origin of this fact. Since it can be already clearly seen in 
perturbation theory, we consider the perturbative Einstein-Hilbert flow f|44l) . The linearized 
flow is given by the matrix 

M= \ dK dd = 

\ OA 8G / 

At the Gaussian FP this matrix becomes 

( -2 L4i(0) 
I 2 
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-2 + BxG 



, XfidM 
r<2dBi_ 




(72) 




Cutoff type la g Cutoff type lb 




-0.2 -0.1 0.1 0.2 0.3 0.4 0.5 -0.2 -0.1 0.1 0.2 0.3 0.4 0.5 



FIG. 7: The flow near the perturbative region with cutoff's of type la and lb. The boundary of 
the shaded region is a singularity of the beta functions. The curves in light color are "classical" 
trajectories with constant AG. 



Cutoff type II 




FIG. 8: The flow near the perturbative region with cutoffs of type II and III. The boundary of the 
shaded region is a singularity of the beta functions. 

which has the canonical dimensions of A and G on the diagonal, as expected. However, the 
eigenvectors do not point along the A and G axes. At the Gaussian FP the "attractive" 
eigenvector is in the direction (1,0) but the 'repulsive" one is in the direction (Ai(0)/4, 1). 
The slant is proportional to ^4i(0) and can therefore be seen as a direct consequence of 
the running of the vacuum energy. This fact has a direct physical consequence: it is not 
consistent to study the ultraviolet limit of gravity neglecting the cosmo logical constant. One 
can set A = at some energy scale, but if G ^ 0, as soon as one moves away from that 
scale the renormalization group will generate a nontrivial cosmological constant. This fact 
persists when one considers the renormalization group improved flow. 

Let us now come to the nontrivial FP. We begin by making for a moment the drastic 
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approximation of treating A\ and B\ as constants, independent of A (this is the leading 
term in a series expansion in A). Thus we consider again the perturbative Einstein-Hilbert 
flow (144j) . In this approximation the flow can be solved exactly: 




(2Ap- |Aigp(l -e 4t ))e- 
2 + 5iG (l -e 2t ) 



The FP would occur at A* = — Ai/4Bi, = —2/Bi, at which point the matrix (|72|) 
becomes 



It has real critical exponents 2 and 4, equal to the canonical dimensions of the constants 
= 2ZA and = —Z. This should not come as a surprise, since the linearized flow 



dimensions, and the eigenvalues are invariant under regular coordinate transformations in 
the space of the couplings. So we see that a nontrivial UV-attractive FP in the A-G plane 
appears already at the lowest level of perturbation theory. It has the form shown in figure 



All the differences between the perturbative Einstein-Hilbert flow and the exact flow are 
due to the dependence of the constants A\ and B\ on A, and in more accurate treatments to 
the RG improvements incorporated in the flow through the functions A2, B2, A3, B%. Such 
improvements are responsible for the nonpolynomial form of the beta functions. In all these 
calculations the critical exponents at the nontrivial FP always turn out to be a complex 
conjugate pair, giving rise to a spiralling flow. The real part of these critical exponents 
is positive, corresponding to eigenvalues of the linearized flow matrix with negative real 
part. Therefore, the nontrivial FP is always UV-attractive in the A-G plane. Conversely, 
an infinitesimal perturbation away from the FP will give rise to a renormalization group 
trajectory that flows towards lower energy scales away from the nontrivial FP. Among these 
trajectories there is a unique one that connects the nontrivial FP in the ultraviolet to the 
GauBian FP in the infrared. This is called the "separatrix" . 

One noteworthy aspect of the flow equations in the Einstein-Hilbert truncation is the 
existence of a singularity of the beta functions. In section IV E, when we neglected the 




(75) 



matrix for the couplings and g^ is diagonal, with eigenvalues equal to their canonical 



m 
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Scheme 


A* 









la 


0.1932 


0.7073 


0.1367 


1.475±3.043i 


la - 1 loop 


0.1213 


1.1718 


0.1421 


1.868±1.398i 


lb with fr 


0.1715 


0.7012 


0.1203 


1.689± 2.486i 


lb with fr - 1 loop 


0.1012 


1.1209 


0.1134 


1.903± 1.099 i 


lb without fr 


0.2329 


0.5634 


0.1312 


2.205 ±3.214 i 


lb without fr - 1 loop 


0.2302 


0.7450 


0.1715 


2.430±2.383 i 


II 


0.0924 


0.5557 


0.0513 


2.425±1.270i 


II - 1 loop 


0.0467 


0.7745 


0.0362 


2.310± 0.382 i 


III 


0.2742 


0.3321 


0.0910 


1.752±2.069 i 


III - 1 loop 


0.0840 


0.7484 


0.0628 


1.695± 0.504 i 



TABLE II: The nontrivial fixed point for Einstein's theory in d = 4 with cosmological constant. 

cosmo logical constant, they appeared at some value G c > G*. Now, looking at equations 
(145f52f57f66p . we see that there are always choices of A and G for which the denominators 
vanish. The singularities are the boundaries of the shaded regions in figures [7] and [SJ Of 
course the flow exists also beyond these singularities but those points cannot be joined 
continuously to the flow in the perturbative region near the Gaufiian FP, which we know 
to be a good description of low energy gravity. When the trajectories emanating from the 
nontrivial FP approach these singularities, they reach it at finite values of t and the flow 
cannot be extended to t — > — oo. The presence of these singularities can be interpreted as a 
failure of the Einstein-Hilbert truncation to capture all features of infrared physics and it is 
believed that they will be avoided by considering a more complete truncation. Let us note 
that for cutoffs of type I and II the singularities pass through the point A = 1/2, G = 0. 
Thus, there are no regular trajectories emanating from the nontrivial FP and reaching the 
region A > 1/2. However, for type III cutoffs the shaded region is not attached to the A 
axis and there are trajectories that avoid it, reaching smoothly the region A > 1/2. 

In table [III we collect the main features of the UV-attractive FP for the Einstein-Hilbert 
truncation with cosmological constant for the different cutoff schemes. 
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V. ULTRAVIOLET DIVERGENCES 



The Einstein-Hilbert truncation does not give a closed set of flow equations, in the sense 
that the beta functions of the higher couplings, which have been neglected in the previous 
section, are not zero. So, if we assume that the higher couplings vanish at some initial 
scale, they will immediately appear as one integrates the flow equations. Before discussing 
truncations that involve higher derivative terms, it will be instructive to see, using the 
ERGE, how such terms are generated in Einstein's theory and how this is related to the 
issue of ultraviolet divergences in perturbation theory. 

A cautionary remark is in order here. In perturbation theory, the divergences appear in 
the formulae relating bare and renormalized couplings. We recall that in our approach we 
never talk of the bare action; instead, we follow the flow of the renormalized action Tk as 
k — > oo. In this limit divergences can appear. However, the limit of r& for k — »• oo cannot be 
simply identified with the bare action. Exploring the relation between these two functionals 
would require introducing an ultraviolet cutoff. We are not going to do this here. In the 
following we will simply compare the divergences of to the perturbative ones. 

In the perturbative approach to quantum gravity, the analysis of ultraviolet divergences 
plays a central role. This issue is not so prominent in the modern literature on asymp- 
totic safety, but this does not mean that divergences do not occur. In an asymptotically 
safe theory, the asymptotic behavior of every quantity is dictated simply by dimensional 
analysis. The dimensionless "couplings in cutoff units" g\ defined in (CQ) tend to constant 
values, so the dimensionful couplings ^ must run like k d \ The couplings with positive mass 
dimension diverge, and those with negative mass dimension go to zero. So, for example, 
near the nontrivial FP in the Einstein-Hilbert truncation discussed in the previous section, 
the graviton wave function renormalization Z = (167rG) _1 diverges quadratically and the 
vacuum energy 2AZ diverges quartically at the FP. 

This matches the powerlike divergences that one encounters in perturbation theory when 
one uses an ultraviolet cutoff. However, in the Wilsonian context these divergences have 
a different physical meaning. The parameter k has not been introduced in the functional 
integral as an UV regulator, rather as an IR cutoff, and in any physical application k corre- 
sponds to some externally prescribed scale. So, in the Wilsonian approach the divergences 
would seem to acquire almost a physical character: they are a manifestation of the depen- 
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dence of the couplings on an external parameter, and it should not be too surprising that if 
an input parameter is allowed to tend to infinity also some output could tend to infinity. 

At a deeper level, however, one should take into account the fact that a dimensionful 
quantity does not have an intrinsic value and therefore cannot be observable. In order to 
give a value to a dimensionful quantity q, one has to specify a unit u, and the result of any 
measurement gives only a value for the dimensionless ratio q/u. Near a FP, the divergence of 
a coupling gi with positive mass dimension for k — > oo is just a restatement of the fact that 
gi, measured in units of k, tends to a constant. If we choose another unit u, since u is also 
ultimately expressible in terms of other couplings, it will also be subject to RG flow. Then, 
the limit q(k)/u(k) may tend asymptotically to zero, to a finite limit or to infinity depending 
on the behavior of u. This highlights that the divergence of a dimensionful coupling cannot 
have a direct physical meaning. The only intrinsic (unit- independent) statement that one 
can make about a dimensionful quantity is whether it is zero, positive or negative. 

Furthermore, only dimensionless functions of the couplings have a chance of being ob- 
servable. It is only for such combinations that the theory is required to give unambiguous 
answers, i.e. it is only such combinations that one could expect to be scheme-independent. 
Now, very often the choice of k which is appropriate to a specific experimental setup is not 
entirely unambiguous. Rather, k sets a characteristic scale of the problem and is usually 



known only up to a factor of order one (see for example the discussion of equation fl6l) or 16] 
for some concrete examples in a gravitational context). The reason why k is nevertheless a 
useful quantity in practice is that dimensionless functions of k tend to depend weakly on k 
and so an uncertainty of order one in the value of k produces only a very small uncertainty 
in the value of the observable (think for example of the logarithmic running of gauge cou- 
pling constants). However, if a coupling is dimensionful, the corresponding dimensionless 
variable cji depends strongly on k, so one should not always expect the value of cji at a given 
scale to be precisely defined. In particular, one should not expect the value of cji at the FP 
to be scheme-independent. 

These expectations are confirmed in the previous treatment of gravity in the Einstein- 
Hilbert truncation. In two dimensions Newton's constant is dimensionless; its value at the 
nontrivial FP (namely zero) and the slope of the beta function are scheme-independent. In 
four dimensions the position of the nontrivial FP in the A-G plane is scheme-dependent. 
However, for all cutoff schemes that have been tried so far A and G are always positive: the 
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existence of the FP and the sign of the couplings are robust features of the theory. In fact, 
in order for A* and G* to be zero, one would have to find a cutoff function such that the 
Q-functionals in front of the heat kernel coefficients B 2 and B are zero. Since and d t Rk 
are positive functions, one sees from ()A14j) that no such choice exists. 

The dimensionless combination AG is related to the on-shell effective action and is known 
to be gauge-independent f|. Numerical studies have also shown that the value of AG at the 
FP is only very weakly dependent on the cutoff function, much less so than the values of A 
and G. It is expected that this residual weak dependence is only an effect of the truncation. 
In fact, it has been argued in 10] that the weakness of this dependence is a sign that the 
Einstein-Hilbert truncation must be stable against the inclusion of further terms in the 
truncation. We will see in section VII that this is indeed the case. 

The UV behavior of dimensionless couplings (i.e. those that are marginal in power count- 
ing) requires some additional clarification. According to the preceding discussion, they have 
a chance of being physically measurable and the existence of a FP requires that they have a 
finite limit. On the other hand in perturbation theory they generically present logarithmic 
divergences. How can these two behaviors be reconciled? It is necessary here to distinguish 
two possibilities: the limit could be finite and nonzero, or it could be zero. If in a certain 
theory all couplings have the former behavior, then there cannot be any logarithmic diver- 
gences. On the other hand if the coefficient g of some operator diverges logarithmically, 
its inverse will go to zero. So if the coupling is the inverse of g, it is asymptotically free. 
This is what happens in Yang-Mills theories, where the (square of the) asymptotically free 
Yang-Mills coupling is the inverse of the coefficient of F 2 . 

In the derivative expansion of four dimensional gravity it is the terms with four derivatives 
of the metric that have dimensionless coefficients. We can parametrize this part of the action 
as follows: 



\c 2 + -R 2 + -E + -V 2 R 
2A 4 p r 



(76) 



where we use the notation introduced in (121j) . The question then arises, what is the asymp- 
totic behavior of these couplings, in particular what is the behavior of A and £? This issue 
can be addressed at various levels, the most basic one being: if we start from Einstein's 
theory, do we encounter divergences proportional to these terms? 

It was shown early on by 't Hooft and Veltman 33] using dimensional regularization that 
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(neglecting total derivatives) the one loop effective action contains the following simple pole 
divergence 



—R UV W V + —R 2 
20 120 



(77) 



Can this result be seen within the ERGE? Let us return to the Einstein-Hilbert truncation. 
In the previous section we have expanded the r.h.s. of the ERGE using the heat kernel 
formula (1A10I) and retained only the first two terms, which are sufficient to give the beta 
functions of the cosmological constant and Newton's contant. Keeping the same inverse 
propagators, we can now consider the next terms in the heat kernel expansion, which will 
give the beta functions of A, £, p, r or more precisely the dependence of these beta functions 
on Newton's constant and on the cosmological constant. We begin by considering a type II 
cutoff; the terms 0(R 2 ) which were not computed in ( 1561) are 
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d t Rk + yRk 
P k - 2A 
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The 64 coefficients for the relevant operators can be computed using equation (]A8j) and the 
traces in (1321): 
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In order to compare to the 't Hooft-Veltman calculation we use the one loop approximation 
to the ERGE, which as explained in section II consists in neglecting r/, and we also set A = 
in (1781 . This gives a contribution to the ERGE equal to 



dt 
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(81) 



where in the last step we have used the identity C 2 = E + 2 (R fll/ R fJjl/ — |-R 2 )- From here one 
can directly read off the beta functions. When one then solves for the flow, in the limit of 
large k each of these couplings diverges logarithmically, with a coefficient that can be read 
off (JHTl) . Recalling that for 't Hooft and Veltman ^ corresponds to ^logA^y, where A uv 
is an UV cutoff, we find agreement with their result in the topologically trivial case. On the 
other hand, if we assume R^ v = 0, which is the on shell condition in perturbation theory 
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without cosmological constant, (1ST]) agrees with the one loop divergence computed in [34]. 
This provides an independent check on the coefficient of the Euler term. Furthermore, our 
calculation shows that all these terms, being proportional to Q ( ) , are independent 



k 

of the choice of the profile function Rk. We have also verified that calculating these terms 
with cutoff's of type la and III leads to the same results. Thus, these divergences are indeed 
independent of the cutoff scheme 9 . 

In the case of cutoffs of type lb there is a subtlety that needs some clarification. With 
these cutoffs (whether one performs a field redefinition, as in section IV B, or not, as in 
Appendix D) it is only possible to perform the calculation on Euclidean de Sitter space (a 
4-sphere). This provides a check on a single combination of the terms appearing in (1761) . 
Specializing to the sphere and using that the volume of the sphere is 3847r 2 /i? 2 , ( JHTj) becomes 



dt 



419 

= if ■ < 82 > 

^ 2 



Using type lb cutoffs one has to pay special attention to the contribution of some of the 
lowest modes in the traces. This is due to the fact that the traces over vector and scalar 
modes may have a prime or a double prime, meaning that some modes have to be left out. 
When evaluating the ERGE without redefining the fields £ M and a, as in equation (Dl), the 
isolated modes give an overall contribution -14, which adds up to the contribution of the 
rest of the spectrum, which is equal to ^^p, to give the correct result. This is an important 
consistency check on the expression (Dl). A similar result holds for the calculation when 
the fields ^ and a are redefined, as in equation (1511 . 

It is interesting to consider also the case when A ^ 0. In the case of a type II cutoff, it 
appears from flTHj) . with rj = 0, and using (1A22IIA23IIA24I) . that the logarithmic divergence 
will be the same as in the case A = 0. This is due to the fact that expanding the fraction in A, 
terms containing A give rise to power-like divergences, so only the leading, A-independent 
term contributes to the logarithmic divergence. The same will be the case for type I cutoffs, 
since again A only appears in denominators. On the other hand for a type III cutoff, ( 1781) 
should be replaced by 



\Qo (?jr) tr&4(A 2 - 2A1) - Q ) tMA (gh)/ 



(83) 



In these calculations wc stick to the dc Donder gauge with a = 1. 
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Then, assuming R^ v = Ag^, which is the on shell condition in perturbation theory with 
cosmological constant, 



dt 



■^R 2 



On the 4-sphere this gives, instead of (1821) 

dT k 



dt 



—R R^ pa - —A 2 
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(84) 



^R 2 



15 



(85) 



These results agree with those obtained in 



35]. Note therefore that the A-dependent con- 



tributions to the logarithmic divergences are scheme- dependent. This should not come as a 
surprise, in view of the discussion above. 

We observe here for future reference that certain authors subtract from the ghost term the 
contribution of the ten lowest eigenvalues of —V 2 on vectors, which correspond to the ten 



Killing vectors of S 4 
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37 



see also 



39 



40J for a discussion. This amounts to putting a 



prime on the ghost determinant (e.g. putting an extra prime on the last two terms in (1501 ). 
This can be motivated by the observation that Killing vectors generate global symmetries, 
and global symmetries are not gauge transformations. Since the ghost contribution has a 
minus sign, with a cutoff of type II this corresponds to adding to ( 1851) the term lO^p^ 
evaluated on the lowest eigenvalue of the operator —V 4 — j, which is equal to zero. This is 



equal to 10 



9tR k {0) 
Pk(0) 



20. Thus, (1851) would be replaced by 
dT k 



dt 



-R? 



329 
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(86) 



We will return to this issue in section VII, when we define the ERGE for f(R) -gravity. 

Having found the expected agreement with earlier one loop calculations based on the 
Einstein-Hilbert action, the next level of sophistication would be to include the terms in 
(ITS"]) in the truncation, as is required in a more accurate approximation to the exact flow. 
We will discuss this in the next section. In the rest of this section we shall discuss the 
possible appearance of divergences that are cubic or of higher order in curvature, keeping 
the kinetic operator that comes from the Einstein-Hilbert action. 

Among higher powers of curvature, of particular interest is the term cubic in the Riemann 
tensor. In perturbation theory, the divergence (ITT]) can be absorbed into a redefinition of 
the metric and therefore does not affect the S matrix: pure Einstein theory is one-loop 
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renormalizable. The first divergence in the effective action that cannot be eliminated by a 
field redefinition in perturbation theory is proportional to R tiU pa R pcr al3 . The coefficient 



of this term was calculated in 



41] at two loops. Can this divergence be seen in the Einstein- 



Hilbert truncation of the ERGE, in the same way as we have seen the 't Hooft-Veltman 
divergence? Let g be the coefficient of this operator in the Lagrangian. In the one loop 
approximation, neglecting A and using a type II cutoff, the beta function of g will be 
proportional to Q-\ = k~ 2 Q where Q is a scheme-dependent dimensionless number. 

As explained in Appendix A, one can choose the cutoff scheme in such a way that Q = 
(this is the case, for example, with the optimized cutoff used in this paper). Then, the 
coupling g will have a FP at g* = 10 . These facts are not surprising: they are a reflection 
of the fact that the ERGE has the structure of a one loop RG equation, and of the absence 
of a Riemann-cube divergence at one loop in perturbation theory. The Goroff-Sagnotti 
counterterm can only be seen in perturbation theory at two loops. On the other hand, if we 
truncate the ERGE at higher order, for example including all terms cubic in curvature, it is 
expected that the beta function of g, though scheme-dependent, cannot be set to zero just 
by a choice of cutoff. There are then two possibilities. The arguments given in the end of 
section III suggest that a FP will exist for all terms in the derivative expansion, including 
the Riemann-cube term; in this case the Goroff-Sagnotti divergence would be an artifact 
of perturbation theory. Alternatively, it is possible that the FP will cease to exist when the 
Riemann-cube term (or some other high order term) is added to the truncation. It is not 
yet known whether this is the case or not; in the conclusions we give some reasons why we 
believe in the former alternative. 



VI. CURVATURE SQUARED TRUNCATIONS. 

In the preceding section we have seen that starting from the Einstein-Hilbert action, 
the RG flow will generate terms quadratic in curvature. In particular, we have discussed 
the way in which these terms diverge logarithmically when k — > oo. It is therefore not 
consistent to neglect these terms, and a more accurate treatment will take into account the 
contribution of these terms to the r.h.s. of the ERGE. In other words, we should include the 

10 if we chose another cutoff such that Q ^ 0, g will go asymptotically to zero; however, as discussed before, 
this is a unit-dependent statement and does not have much physical meaning. 
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terms (176|) in the truncation. The resulting theory is perturbatively renormalizable 58| but 
has problems with unitarity. Furthermore, in the perturbative treatment the cosmological 
constant is set to zero, something we know can only be imposed at a given scale. The 
corresponding Wilsonian calculation has been done within the one loop approximation, and 
has been briefly reported in [ljj. Here we review and extend those results. 

We take the action Tk as the sum of the Einstein-Hilbert action (1281) and the terms 
7G\i . The linearized wave operator is now a complicated fourth order operator. In order to 
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43 



44 



451 ] it is convenient to choose a gauge fixing of the 



Sgf = / d*x^ X ,Y^ Xl 



simplify its form, following 
form 



where \v = V^fym, + /3V,yh (all covariant derivatives are with respect to the background 
metric) and 

a 

The ghost action contains the term 



where 



S c = / d*xJgC v {A gh )\C» 



(A gh )\ = -Sy 2 - (1 + 2(3) \/y + R\ 



as well CIS db "third ghost" term 



Sh = \ I 



due to the fact that the gauge averaging operator Y depends nontrivially on the metric. We 
follow earlier authors 44] in choosing the gauge fixing parameters a, (3, 7 and 5 in such a 
way that the quadratic part of the action is: 



where 



\ J d A x^ 5gKA^5g 



A (4) = (-V 2 ) 2 + V A V P V A + U 



(87) 



(88) 



441 ] whose notation we mostly 



For details of the operators K, V and U we refer the reader to 
follow. We choose type III cutoffs: 

R|(A( 4 )) = KR^(A^) ; Rt(A {9h) r u = 5^(A {gh) ) ; ^(Y)^ = g^Rf{Y) . 
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For higher derivative operators we will use a generalized optimized cutoff of the form 



R 



(ak n — z)8(ak n — z), with a = 1 unless otherwise stated. 



We restrict ourselves to the one-loop approximation, as explained in previous sections. 
Furthermore, only the contributions of the heat kernel coefficients up to B4 will be taken into 
account. In this way, we will essentially neglect all RG improvements. The beta functions 
of the dimensionless couplings appearing in (1761) turn out to be: 
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They form a closed system and agree with those calculated in dimensional regularization 
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44 



451 ] . It is convenient to define new couplings ui and 9 by £ = —3X/u and p = X/8. 



In this way 1/A is the overall strength of the terms quadratic in curvatures, and 9 and u 
give the relative strength of the different invariants. Then, the beta functions become 

1 25 + 1098 u + 200 a; 2 



U 
0e 



(4tt) 2 
1 7(56 
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(4tt) 
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(90) 



The coupling A has the usual logarithmic approach to asymptotic freedom, while the other 
two couplings have the FP values a;* m (—5.467, —0.0228) and 6* ~ 0.327. Of the two roots 



for uj, the first turns out to be UV-repulsive, so the second has to be chosen 
The beta functions of A and G are: 
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45|. 
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256ttGu 2 
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86lu + 40a/ 
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aA 2 A + q(u)GA 

b4ir z uj it 
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(91) 



(4vr) 2 12a; 

: (83 + 70a; + 8a; 2 )/187r. The first two terms in each beta function exactly 
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45], the remaining ones are new. The origin of the new 



reproduce the results of 43 
terms will be discussed below. 

To study the flow of A and G, we set the remaining variables to their FP values u 



00* 



6 = 6*, and A 



A* = 0. Then, the flow takes the form of the perturbative Einstein-Hilbert 
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flow g3D, with A x = 4/ V, B 1 = -g* = « -1.440, A 2 = 5 2 = (see figured]). It has 

two FPs: the Gaussian FP at A = G = and another one at 

1 ~ 2 

A* = « 0.221 , = — ~ 1.389 . 

Trq* q* 

Like all one loop flows with constant A\ and B\, the critical exponents are 2 and 4. This, 
however, is due to the approximation. If we could take into account the contribution of the 
heat kernel coefficients B$, B$ etc., they would contribute terms of order A 2 and higher to 
the beta functions. We expect that these terms would produce a complex conjugate pair of 
critical exponents, and the corresponding spiralling flow. In this connection see also 55]. 

In the preceding calculation we have used, besides the truncation to four-derivative terms, 
also the one loop approximation, and contributions coming from the heat kernel coefficients 
Bq and higher have been neglected. These are also the approximations made in earlier 



perturbative calculations 42, Q, |44j, |45j, so it is instructive to understand the origin of 



the additional terms in ( 191~|) . which are essential in generating the nontrivial FP. The beta 
functions were originally derived as coefficients of 1/e poles in dimensional regularization, 
which correspond to logarithmic divergences in the effective action. In a heat kernel deriva- 
tion these terms are given by the B^ coefficient. In the old calculations only these terms 
were retained. The new terms that we find come from the B^ and Bq coefficients, which in 
a conventional calculation of the effective action would correspond to quadratic and quar- 
tic divergences. Such terms are discarded in dimensional regularization, but we see that 
proceeding in this way one would throw away essential physical information. In order to 
keep track of this information in dimensional regularization one would have to take into 
account the contribution of the pole occurring in dimension 2. This can be done with the 
2 + e expansion, and we have already discussed the way in which the FP then appears. The 
situation is entirely analogous to the Wilson-Fisher FP in three dimensional scalar theory, 
which can be seen either using a cutoff regularization or, if dimensional regularization is 



used, in the 4 — e expansion [46| . It is appropriate to stress once more that our ' Wilsonian" 
calculation of the beta functions does not require any UV regularization. Accordingly, there 
are no regularization/renormalization ambiguities; the only ambiguity is in the choice of the 
cutoff procedure, but we have already seen that no choice of Rk could remove the B2 and 
Bq terms. 

Having discussed the beta functions in the higher derivative truncation of pure gravity, 
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it is a simple exercise to add to them the contributions of minimally coupled matter fields, 
which were discussed in section III. From (189]) and fl2TT) we find 

R - 1 133 \2 o\2_(4) 

& = TZ^^ ~ , (92) 
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2880vr 2 2 



= ( — uq rip — 31um) . (93) 

p 2880tt 2 V 2 2 J K J 

For these couplings, the new terms simply change the direction and speed of the logarith- 
mic approach to asymptotic freedom. In particular, the ratios of the couplings approach 
asymptotically the following FP values 

U * = 250 (~ 549 ~ 9607r2a A 4) ±\/9216007r 4 (4 4) ) 2 + 1054080vr 2 4 4) - 5760004 4) vr 2 + 296401^ 



.2^(4) 

6* 



8 49 + 180vr 2 a), 



9 133 + 3207r 2 ai 4) 

On the other hand, the beta functions of A and G are modified by the addition of Sna^G + 
^-na^'AG and 32-na^G 2 respectively, where, using a type II optimized cutoff, 



1 



a (0) = ^ (n s - 4n D + 2n M ) 
1 



a 



(2) 



(n s -2n D -4n M ) . (94) 



96vr 2 

Let us consider again the flow in the A-G plane. Putting uo and 9 to their FP values and 
taking the limit A — > 0, the flow takes again the form (|43p . with the coefficients A\ = 
4/vr + 16vra(°\ B 1 = -q{uj*) + 16W 2 ), A 2 = B 2 = 0. 

This calculation is of some interest for the following reason. Recall that whereas applying 
perturbation theory to pure Einstein theory, terms proportional to RapR pl/ and R 2 can be 



absorbed by a field redefinition, this is no longer so when matter fields are present 33]. 



Thus, gravity coupled to matter is nonrenormalizable already at one loop. One may fear 
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that the FP ceases to exist when one includes in the truncation terms that correspond 
to nonrenormalizable divergences in the perturbative treatment of Einstein's theory. The 
preceding calculation shows, at least at one loop, that this is not the case. We will further 
comment on the significance of this point in the conclusions. 



VII. TRUNCATION TO POLYNOMIALS IN R 



Following the scheme outlined in the previous two sections, we now extend the formal- 
ism to handle truncations in the form of the so-called "/(-R)~gravity" theories, where the 
Lagrangian density in (JSj) is a function of the Ricci scalar only. Such theories rose strong 



interest recently in cosmological applications 43]. At one loop, the quantization of these 



theories has been discussed in [48[. Here we analyze the RG flow of this type of theories 
where / is a polynomial of order n < 8. 

The (Euclidean) action is approximated by 

n „ 

r*[$] = 9^ k ) / ddx V9 Rt + 3 gf + S gh , (95) 

i=0 

where $ = {h ,,,,, C /t ,, C„} and the last two terms correspond to the gauge fixing and the 



ghost sector 



18 



321 ] . The gauge fixing will have the general form 

S gf = \J d d x^X^ X u (96) 

where Xv — V^/i^— ^ £ V I ,/i'i and G^ v = (a + (3V 2 ) g^. In the following four subsections we 
give the second variation of the action and the ghost action, then we insert these expressions 
into the ERGE and finally we discuss the choice of gauge. 



A. Truncation Ansatz 

The second variation of J & d x^/gf (R) gives 
S 2 j d d x (y/gf (R)) = J d d x y/gl \f" (R) (^V a V^WV - 2W 2 V a V p h. 



a/3 

+h(V 2 ) 2 h - 2V a V p h aP R^W v + 2V 2 hR flu h^ + R^R aP h^h afi ) 
1 
2 



+\f (R) (h a(3 v 2 h al3 + hsr^vhr - 4h a/3 v a whP + 2h a pWV a h? 



-hV 2 h + AR^h u „ - 2R ilv W v h) + / (R) {-\K V W U + 
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All covariant derivatives are with respect to the background metric, the trace is abbre- 
viated as h, prime denotes derivative with respect to R. It is already assumed here that the 
curvature tensor of the background metric is covariantly constant. To this one has to add 
the gauge fixing terms (IDT?]) . In order to diagonalize the complete expression, we choose a 
(Euclidean) de Sitter background and decompose into tensor, vector and scalar parts as we 
did in section IVB for the cutoff of type lb. Then one obtains for the tensor part 
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and for the scalar part (which contains a nontrivial mixing between h and a) 
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We have dropped the subscript k for typographical clarity. 
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B. Ghost terms 



The Fadeev-Popov ghost consists of two parts. It is calculated in the usual way from the 
variation of the gauge condition and the generators of gauge transformations by 



(100) 
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From the infinitesimal gauge transformation of the gravitational field, 5h^ v = V^e,, + V^e^, 
the variation of the gauge condition under gauge transformations is given by 

a 

= V\ + + d ~ 2 ~ 2p V u V^. (101) 



This gives 



S c = J d d x^C^{a + (3V 2 ) 



^V 2 + R\ + {d ^ 2p V v, 
a 



C u (102) 



where C M and C 1 are the ghost and anti-ghost fields. 

As we want to treat also higher-derivative gravity, it is natural to assume the operator 
can contain derivatives (for (3^0). In that case, one obtains a nonconstant square root 
of a determinant in the Fadeev-Popov procedure which on exponentiation gives rise to the 



so-called third ghost term 



45) 



S b = ~ J tfxy/jjWG^W. (103) 

For (3 = this gives just a constant factor which can be absorbed in the normalization of 
the functional integral. The full ghost action is then S g h = S c + St,. The ghost, anti-ghost, 
and third ghost fields are also decomposed into transverse and longitudinal parts as in (1471) . 
The operators acting on these fields are: 



(2) _ _2(d-l-p) / , fl /_ 2 fl\ w_ 2j D \_ 2 



= (« + /3V 2 )<r 

r£ } = -(a + ^(v 2 + ^)v 2 . (104) 
Finally, the decomposition of the ghosts gives rise to Jacobians involving the operators 

J c = -V 2 ; J b = -V 2 . (105) 



C. Inserting into the ERGE 

We choose a cutoff of type lb. The inverse propagators (I97II98II99|I104I) are all functions 
of —V 2 . Then, for each type of tensor components, the (generally matrix- valued) cutoff 
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function Rk is chosen to be a function of —V 2 such that 

r (2)(_v 2 ) + R k (-V 2 ) = T^(P k ) , (106) 



where Pk is defined as in ( TTTj) for some profile function Rk. Inserting everything into the 
ERGE © gives: 

dTk _ 1 Tr Tt R hT h T 1 rp / Tt R tt 
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__Tr' dt^Jv _ ^rp r // dt^Js 
2 Wjy + R Jv 2 {0) Js + Rj s 



-Tr' (0) dt Jc - -TiLs dt Jb . (107) 



dp i d p 

+ i? Jc 2 (0) J b + i 

The first three lines contain the contribution from the metric fluctuation h^, which has been 
decomposed into its irreducible parts according to (|46|) . Note that the trace over the scalar 
components is doubly primed, since the first two modes of the a field do not contribute to h^ v . 
However, the first two modes of h do contribute, and their contribution to the trace is added 
separately in the third line. The fourth and fifth lines contain the contributions of the ghosts 
and the third ghost, each decomposed into transverse and longitudinal parts. Note that the 
first mode of the scalar (longitudinal) parts is omitted, as it does not contribute to C M and 
respectively. The sixth line is the contribution of the Jacobians of the transformation 
(146]) . These have to carry the same number of primes as the fields in ( T46]) . The same is valid 
for the traces over the Jacobians resulting from the split into vector and scalar parts of the 
ghost and third ghost field, given in the last line. Eliminating the Jacobians by a further 
field redefinition, as in section IVB, would produce some technically undesirable poles in 
the heat kernel expansion. For this reason we shall proceed as in Appendix D and explicitly 
retain the Jacobians. 
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In equation (j!07p we have used the convention for the primes that was used in section IV 
B and in Appendix D. This differs from the calculation we did in [ijj] in having less primes 
in the traces over the ghosts. We will discuss this point in some detail in the next section. 



D. Discussion of gauge choices 

In the gauge fixing term (1961) we have allowed in principle three gauge fixing parameters 
p, j3, and a. The choice p — 1, (3 — corresponds to the familiar de Donder gauge; in 
the discussion of the Einstein-Hilbert truncation we have further chosen the gauge fixing 
parameter a = Z, which produces a minimal kinetic operator. Other values of a have been 
treated in [h]]. To avoid the issue of the RG running of a, the limit a/Z — > oo is sometimes 
invoked. A gauge fixing with (3 ^ contains terms with four derivatives and is natural in 
higher derivative gravity. We will only discuss gauge choices where either a or (3 are nonzero, 
and not both simultaneously. We will call these "a gauges" and u /3 gauges" respectively. 
Note that p and (3 are dimensionless but a has dimension of mass squared. There are two 
ways of turning it into a dimensionless parameter: the first is to proceed as with all other 
couplings and define at = ak~ 2 ; then we can set a = 1/a. The second is to proceed as in 
(]29p for the Einstein-Hilbert truncation and set a = Z/a. In the following, when we use a 
gauges we will always adopt the first method; the second method yields similar results (up 
to a rescaling of a) and will not be reported. We will also use the definition b = 1/(3. We 
will always neglect the RG running of the dimensionless gauge parameters p, a and b. 

To reach the highest degree of simplification, a convenient gauge choice is to set p = 
and then take either (3 = and a — > oo or a = and (3 — > oo. We will now show explicitly 
how this simplification works in the (3 gauges, treating tensor, vector and scalar components 
separately. To write the formulae in a more compact form, in this section we will denote 
A = — V 2 and we define the following shorthands 



The transverse traceless tensor part is gauge independent and therefore is not affected by 



AW = A - Pi 




the gauge choice. The vector part receives contributions from £, c 



;T C Tfl b hr- 
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1 dt [-g (P k - A)] 1 , dt (P k - A) 
+ 2 (1) (a-W) " 2 (1) pf ■ 
In the first term, looking at (|98|) one sees that in the limit (3 — > oo only the gauge fixing 
term matters. In this limit the first term becomes simply 

5 (1> ^fT~s w ~* 



Treating the other three terms in the same way, several simplifications occur. However, 
one has to pay attention to the fact that some of the traces are primed and some are not. 
Therefore, in the simplifications the contributions of some isolated modes are left out. Using 
the optimized cutoff and specializing to d = 4, the final result for the vector terms is 

where the argument of the last two terms, j, is the first eigenvalue of the Laplacian on 
transverse vectors (see table WW) . 

Let us now come to the scalar part. It receives contributions from h, a, c, c, b, Js, J c 
and J;,. The contribution of h and a is given by the second line in eq. f l 1 071) . One sees from 
(|99|) that when we set p = 0, only T a <r and R aa contain (3. Therefore in the limit (3 — > oo 



these terms become 

The second term is equal to 

l m „ 9 t i? CT(T 1 .. / 9ii? fc n d t R k d t R k 



The longitudinal ghost fields c and c give a contribution 



t (P fc P fc (3) P fc (4) -AA( 3 >A( 4 >) ( d t R k d t R k d t R k 

(0) ftP^pW " (0) I Pf + Pf ft 



The third ghost and the Jacobians together contribute 



1 , ^(P fc Pf -AAW) _ 1 „ d t (P k P^ - AA (3) ) 

2 W p P (4) 2 (0) P pP) 

<9,(P fc - A) 1 ft(ft-A) 

+ 117 m — — 1*7 



{0) P fc 2 Xi W P fe 




// ( ^i-^fc . <%P fc \ <9*Pfc 1 , d t R k 

2 (0) ^pr + "p7; + (0) ~fT ~ 2 (0) ~a~ 
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Thus the scalar terms together give 

The last two terms are evaluated on the second eigenvalue of the Laplacian on scalars, ^. 
Finally we can collect the tensor, vector and scalar contributions to obtain 

dT k 1 f d t P k f + (P k - A)d t f } 1 , <9 t P fc 1 ., <9 t P fc . . 

2 (2) \ (P k -f)/' + / J 2 (1) P fe -f 2 (0) P fc -f 

1 „ f 5 t P fe (/' + 6(P fe - f )/") + (P k - A)($/' + 3(P fc + A - f P)<9 t f 



4-_Tr" 1 ' K 3 ' ' \ /c — /v-^ i V 1 — 3 

2 | /+(Pfc -|i?) f _ 3m _|)2 



where £ is the contribution of the isolated modes. In the gauge p = 0,a = 0,/3 — > oo it is 

^P fc (|) <9fP fc (f) d t R k (%) 5<9 t P fc (f) , 

^(f) n (4) (f) n (3) (!) 2 ^i 4) (f) 

The calculation in the gauge p = 0, a — > oo, /5 = proceeds in a similar way. The final 
result is the same except for the isolated modes, which give 

d t R k (§) <9 t P fe (#)P<: 6) (#) d t R k (%) , . 

£ = -10 ' - 10 Us; , q fc , Vs; +5 V B 3; . 112 

n (4) (f) ^(f)n (3) (f) P ^(f) 

Using the formulae in the appendix for the trace evaluation, writing R = k~ 2 R and 
/ = k~ 4 f, f = k~ 2 f, f" = f", with the optimized cutoff this equation becomes 

dT k 384vr 2 f 1008(511P 2 - 360P - 1080) 2016(607P 2 - 360P - 2160) 
dt ~ 30240P 2 \ P-3 P-4 

+2Q (31 IP 3 - 126P 2 - 22680P + 45360) (dj' + 2f - 2Rf") - 252(P 2 + 360P - 1080)/' 
+ 3/ -(P-3)/' 

+ [l008(29P 2 + 360P + 1080)/' + 4(185P 3 + 3654P 2 + 22680P + 45360) (d t f + 2f - 2Rf 
-2016(29P 3 + 273P 2 - 3240)/" - 9(181P 4 + 3248P 3 + 15288P 2 - 90720) (dj" - 2Rf"' 

(/" (P - 3) 2 + 2/+ (3 - 2P)f) J + £ . (113) 



where 



10(P 2 -20P + 54)P 2 „ 10(llP-36) 

or £ = = (114) 



P 2 - 7P + 12 P 2 - 7P + 12 
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in the (3- and a-gauge respectively. The (dimensionful) beta functions can be extracted 
from this function by taking derivatives: 

d gi 1 d { 1 dT k 



dt i\dB>V dt 



(115) 



R=0 

where V — J d^x^fg. This we have done using algebraic manipulation software. 

As already mentioned in section V, different authors have different prescriptions for the 
treatment of zero modes of the ghost operator. In [jj]] we have taken the point of view that 
since in the gauge p = 0, a = 0, (3 oo there is a one-to-one correspondence between the 
modes of the unphysical components £ M and a and those of the ghost field, and since the 
ghost contribution is supposed to cancel the contribution of the gauge degrees of freedom 
of the field, the cancellation occurs mode by mode, so that the trace over the vector part of 
the ghost has to have a prime and the trace over the scalar part of the ghost has to have a 
double prime. A similar argument applies to the third ghost. Finally, when one makes this 
choice for the ghosts, then also the Jacobian determinants J c and J& have to have a double 
prime. So, altogether, all vector traces would have a prime and all scalar traces would have 
a double prime. This amounts to putting E = 0. In the next sections we shall begin by 
giving the results with this definition of the traces. Later we shall also describe the effect of 
having £ as in equation (]114p . 



E. Results 

We can now state our results. Table II III gives the position of the nontrivial FP and table 
IIVI gives the critical exponents, for truncations ranging from n — 1 (the Einstein-Hilbert 
truncation) to n = 8. The same information is shown graphically in figure ([9]). 

Some comments are in order. First of all, we see that a FP with the desired properties 
exists for all truncations considered. When a new coupling is added, new unphysical FP's 
tend to appear; this is due to the approximation of / by polynomials. A similar phenomenon 



is known to occur in scalar theory in the local potential approximation |49l . l50|. However, 
among the FP's it has always been possible to find one for which the lower couplings and 
critical exponents have values that are close to those of the previous truncation. That FP 
is then identified as the nontrivial FP in the new truncation. 

Looking at the columns of Tables II I II and [TV] we see that in general the properties of the 
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n 


A* 


G* 


A^G* 


10 3 x 










9o* 


9i* 


92* 


93* 


9a* 


95* 


9&* 


97* 


98* 


1 


0.1297 


0.9878 


0.1282 


5.226 


-20.140 
















2 


0.1294 


1.5633 


0.2022 


3.292 


-12.726 


1.514 














3 


0.1323 


1.0152 


0.1343 


5.184 


-19.596 


0.702 


-9.682 












4 


0.1229 


0.9664 


0.1188 


5.059 


-20.585 


0.270 


-10.967 


-8.646 










5 


0.1235 


0.9686 


0.1196 


5.071 


-20.538 


0.269 


-9.687 


-8.034 


-3.349 








6 


0.1216 


0.9583 


0.1166 


5.051 


-20.760 


0.141 


-10.198 


-9.567 


-3.590 


2.460 






7 


0.1202 


0.9488 


0.1141 


5.042 


-20.969 


0.034 


-9.784 


-10.521 


-6.048 


3.421 


5.905 




8 


0.1221 


0.9589 


0.1171 


5.066 


-20.748 


0.088 


-8.581 


-8.926 


-6.808 


1.165 


6.196 


4.695 



TABLE III: Position of the FP for increasing order n of the truncation. To avoid writing too many 
decimals, the values of have been multiplied by 1000. 



n 


Red! 


Im , &! 


tf 2 


^3 


Rcd 4 




^6 


#7 


^8 


1 


2.382 


2.168 
















2 


1.376 


2.325 


26.862 














3 


2.711 


2.275 


2.068 


-4.231 












4 


2.864 


2.446 


1.546 


-3.911 


-5.216 










5 


2.527 


2.688 


1.783 


-4.359 


-3.761 


-4.880 








6 


2.414 


2.418 


1.500 


-4.106 


-4.418 


-5.975 


-8.583 






7 


2.507 


2.435 


1.239 


-3.967 


-4.568 


-4.931 


-7.572 


-11.076 




8 


2.407 


2.545 


1.398 


-4.167 


-3.519 


-5.153 


-7.464 


-10.242 


-12.298 



TABLE IV: Critical exponents for increasing order n of the truncation. The first two critical 
exponents #o an d "&! are a complex conjugate pair. The critical exponent $4 is real in the truncation 
n = 4 but for n > 5 it becomes complex and we have set $5 = $4. 



FP are remarkably stable under improvement of the truncation. In particular the projection 
of the flow in the A-G plane agrees well with the case n = 1. This confirms the claims made 
in 1(| about the robustness of the Einstein-Hilbert truncation. 

The greatest deviations seem to occur in the row n — 2, and in the columns g% and $2- 
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1 2 3 4 5 




- Rethetal 
theta2 

-theta3 

- Retheta4 
-theta6 
-theta7 

thetae 



FIG. 9: The position of the fixed point (left panel) and the critical exponents (right panel) as 
functions of n, the order of the truncation. 



The value of #2* decreases steadily with the truncation. The critical exponent $2 appears for 
the first time in the truncation n = 2 with a very large value, but it decreases quickly and 
seems to converge around 1.5. This behavior may be related to the fact that g 2 is classically 
a marginal variable. 

The beta function of g 2 due to the Einstein-Hilbert action (in the spirit of section V) was 



considered first in 



51] ; the full truncation n = 2 has been studied in 



our results for the case n = 2 with those of 



When comparing 



52| , one has to keep in mind that they generally 



depend on the shape of the cutoff function. A significant quantity with very weak dependence 
on the cutoff function is the dimensionless product AG. The value 0.12 -j- 0.14 given in 52] 
for AG is very close to the value we find in all truncations except n = 2. Our value for 
g 2 * in the n = 2 truncation has the same sign but is between one half and one third of 
their value, depending on the cutoff function. This is another manifestation of the relatively 
unstable behavior of this variable. The value given in 52| for the critical exponent varies 
in the range 2.2 -j- 3.2 depending on the shape of the cutoff, and is in good agreement with 



52| the critical 



our results, again with the exception of the n = 2 truncation. Finally, in 
exponent $2 has stably large values of the order of 25 with the compact support cutoffs, 
but varies between 28 and 8 with the exponential cutoffs. The values at the high end agree 
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gauge 


A* 


G* 


A^G* 


10 3 x 










9o* 


9i* 


92* 


93* 


94* 


95* 


96* 


97* 


98* 


a 


0.1239 


0.9674 


0.1199 


5.096 


-20.564 


0.153 


-6.726 


-5.722 


-2.981 


1.980 


3.305 


1.631 


(3 


0.1242 


0.9682 


0.1202 


5.103 


-20.548 


0.138 


-6.133 


-4.621 


-1.407 


2.240 


2.207 


0.610 



TABLE V: Position of the FP for n = 8 taking into account the contribution of the isolated modes 
given in (| 1 14j) . To avoid writing too many decimals, the values of have been multiplied by 1000. 



gauge 






tf 2 


#3 


Re$ 4 


Im&4 


??6 


#7 


#8 


a 


2.123 


2.796 


1.589 


-4.212 


-1.107 


5.558 


-7.321 


-9.923 


-12.223 


a 


2.049 


2.511 


1.438 


-3.928 


-0.102 


7.320 


-7.239 


-9.664 


-12.381 



TABLE VI: Critical exponents in the a- and /3-gauge taking into account the contribution of the 
isolated modes given in equation (|114p . 

well with our result in the n = 2 truncation. The shape dependence that is observed with 
exponential cutoffs can be taken as a warning of the truncation-dependence of this quantity. 

Tables |V] and IVII give the position of the FP and the critical exponents in the truncation 
n = 8, using the definition of the traces with less primes (i.e. using for £ the values given 
in equation (11141) ). While the numerical results, especially for some of the higher couplings, 



the overall qualitative picture is not affected. In this connection we also 



do change sensib 

mention that in 53j the same calculation has been independently repeated for n < 6. The 
slight numerical differences between their results and those reported here is due entirely to 
the fact that they define the traces over the Jacobians J c and J& with a single prime; it has 
been checked that when the same definition is used, the results agree perfectly. 



F. UV critical surface 



Possibly the most important result of this calculation is that in all truncations only three 
operators are relevant. One can conclude that in this class of truncations the UV critical 
surface is three-dimensional. Its tangent space at the FP is spanned by the three eigenvectors 
corresponding to the eigenvalues with negative real part. In the parametrization (195]) . it is 
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the three-dimensional subspace in R 9 defined by the equation: 



Q-\ = 

JO 


0.0006 + 0.0682 O + 


0.4635 & + 


O.895O0 2 


9a = 


—0.0092 — 0.8365 g 


- 0.2089^1 + 1.6208^2 


95 = 


-0.0157- 1.2349 O 


- 0.72540! 


+ 1.0175 02 


96 = 


-0.0127-0.6226 g 


- 0.8240 0i 


- 0.6468 02 


97 = 


-0.0008 + 0.8139<7o 


-0.148401 


-2.018102 


98 = 


0.0091 + 1.2543^o + 


0.50850! - 


1.901202 



(116) 

Unfortunately, we cannot yet conclude from this calculation that the operators of 1 ^ with 
n > 6 would be irrelevant if one considered a more general truncation: the beta functions 
that we computed here are really mixtures of the beta functions for various combinations of 
powers of Riemann or Ricci tensors, which, in de Sitter space, are all indistinguishable. 

However, there is a clear trend for the eigenvalues to grow with the power of R. In fact, 
in the best available truncation, the real parts of the critical exponents differ from their 
classical values di by at most 2.1, and there is no tendency for this difference to grow for 
higher powers of R. This is what one expects to find in an asymptotically safe theory 5J]. 

With a finite dimensional critical surface, one can make definite predictions in quantum 
gravity. The real world must correspond to one of the trajectories that emanate from the 
FP, in the direction of a relevant perturbation. Such trajectories lie entirely in the critical 
surface. Thus, at some sufficiently large but finite value of k one can choose arbitrarily 
three couplings, for example 0o, gi, 02 and the remaining six are then determined by (11161) . 
These couplings could then be used to compute the probabilities of physical processes, and 
the relations (11161) . in principle, could be tested by experiments. The linear approximation 
is valid only at very high energies, but it should be possible to numerically solve the flow 
equations and study the critical surface further away from the FP. 



G. Scheme dependence in i? 3 gravity 

The gauge parameters cannot be chosen totally independently. For example, if one takes 
p = 0, a and (3 can be taken to infinity or zero, the limits exist and do commute. If one 
takes instead p ^ 0, either a or (3 has to stay finite. 
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FIG. 10: The position of the nontrivial fixed point (left panel) and the critical exponents (right 
panel) in the R 3 truncation, in the gauge p = and gauge parameter a variable in the range 
— 1024 < a < 64. Note that the scales are logarithmic for a > and a < but not at a = 0. 
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FIG. 11: The position of the nontrivial fixed point (left panel) and the critical exponents (right 
panel) in the R 3 truncation, in the gauge p = 1 and gauge parameter a variable in the range 
0.5 < a < 64. 



It would be very cumbersome to check the stability of the results under changes in the 
gauge fixing parameters in high truncations. However, it is possible to do so for n — 3. As 
we shall see, this is enough to verify that g% is always irrelevant and therefore is a strong 
indication that the dimension of the critical surface is stable under such variations. 

Figures [10] and [12] show the position of the FP and the critical exponents in the gauges 
p = 0, (3 = 0, a variable and p = 0, a = 0, (3 variable, respectively. It appears that one 
can take the limits a — > ±oo and (3 — > ±oo without problems and that the results are 
continuous. 

We have also checked the results in the gauge p — 1, which was discussed extensively 
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FIG. 12: The position of the nontrivial fixed point (left panel) and the critical exponents (right 
panel) in the R 3 truncation, in the gauge p = and gauge parameter b variable in the range 
—32 < b < 256. Note that the scales are logarithmic for b > and b < but not at 6 = 0. 





FIG. 13: The position of the nontrivial fixed point (left panel) and the critical exponents (right 
panel) in the R 3 truncation, in the gauge p = 1 and gauge parameter b variable in the range 
3 < b < 256. 

in the case of the Einstein-Hilbert truncation. Figures [11] and [13] show the position of the 
FP and the critical exponents in the gauges p = 1, (3 = 0, a variable and p = 0, a = 0, (3 
variable, respectively. It appears that when p = 1 one cannot take the limits a — > oo and 
(3 — > oo, so we have limited ourselves to positive values of a and (3. 
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VIII. CONCLUSIONS 



In this paper we have reviewed and extended recent work on the Wilsonian RG approach 
to quantum gravity. The central hypothesis of this approach is the existence of a nontrivial 
FP for gravity, having finitely many UV-attractive directions. Accordingly, most of the 
work has gone towards proving the existence of such a FP. Let us summarize the evidence 
that has been obtained by applying Wilsonian RG methods to gravity (for general reviews 



see also 
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m. 



In order to start from the simplest setting, we have begun by considering the contribution 
of minimally coupled matter fields to the beta functions of the gravitational couplings, 
which can be simply obtained from the heat kernel expansion for Laplace-type operators 
and give beta functions that are just constants (see equations (I24II26P ). Such beta functions 
produce a FP for all gravitational couplings. This result is important for two reasons: the 
first reason is that in the limit when the number of matter fields is very large, this is the 
dominant contribution to the beta functions. Insofar as the number of matter fields in 
the real world is large, and matter couplings can be assumed to be asymptotically free, 
this may be a reasonably good approximation for some purposes. The second reason is 
that the contribution of gravitons to the gravitational beta functions, neglecting the RG 
improvements and considering only the leading term in the expansion in the cosmological 
constant, is essentially of the same form (compare equations (1261) . (1441) and (I9TT) ). One may 
therefore see this as a bare skeleton that can be dressed by taking into account increasingly 
more subtle effects. 

The simplest way of approximating the ERGE consists in truncating the form of the 
effective average action Tk, i.e. retaining only certain operators and neglecting all others. It 
is known from experience with scalar field theory that the most reliable procedure is not to 
truncate on the power of the field but rather to truncate in the number of derivatives. In the 
case of gravity the lowest order of the derivative expansion is the so-called Einstein-Hilbert 
truncation, where one retains only the cosmological constant and Newton's constant. In this 
case one can treat the ERGE without any further approximation. The results are therefore 
"nonperturbative" , in the sense that couplings are not required to be small. In section IV we 
have discussed in detail several ways of implementing the cutoff procedure of the ERGE and 
we have shown that the results are robust against such changes. This adds to earlier studies 
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concerning the dependence of the results under changes of gauge and changes of the cutoff 
profile and increase our confidence that the FP is not an artifact of the truncation. The 
results have been obtained in arbitrary dimensions d. Since in this approach the dimension 
is not used as a regulator, one can follow the position of the FP as a function of d without 
encountering singularities, and compare with the results of the epsilon expansion. We have 
seen that, numerically, the epsilon expansion gives a rather poor approximation at d = 4. 

We have then shown how to recover the perturbative divergences from the ERGE. In 
particular, we have seen that the one loop divergences obtained by 't Hooft and Veltman 
331 ] can be reproduced starting from the Einstein-Hilbert truncation of the ERGE, and that 
they are independent of the cutoff procedure and of the profile of the cutoff functions. We 
have also reproduced the known (scheme-dependent) one loop divergences in the presence 
of a cosmological constant. 

A more accurate treatment of these divergences requires that the terms quadratic in curva- 
ture be retained from the start, i.e. that they are included in the truncations. Unfortunately, 
due to technical complications, it has been impossible so far to treat the curvature-squared 
truncation (which would constitute the second order in the derivative expansion) in the same 
way as the Einstein-Hilbert truncation. The most complete available treatment, which was 
described in section VI, requires that further approximations be made: essentially, one is just 
keeping the lowest terms in the perturbative expansion. Still, the results obtained from the 
ERGE differ from those that had been calculated before using more conventional methods. 
The beta functions of the dimensionless couplings A, £, p (defined in ( 1761) ) coincide with 
those that had been computed previously. In this approximation these couplings are asymp- 
totically free, tending logarithmically to zero from a well-defined direction in (A, £, p)-space. 
However, the beta functions of A and G contain, in addition to the terms that were known 
before, also some new terms that generate a nontrivial FP. When the other couplings are 
set to their FP values, the flow in the (A, G)-plane has the same form as the perturbative 
Einstein-Hilbert flow (Jill . 

In the perturbative approach to Einstein's theory, the one loop divergences are at most 
quadratic in the curvature. But the Euler term is a total derivative and the remaining 
terms vanish on shell, so, if we neglect the cosmological constant, all counterterms can be 
eliminated by a field redefinition, up to terms of higher order. One could suspect that 
in pure gravity the existence of the FP in truncations involving at most terms quadratic 
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in curvature is in some way related to the absence of corresponding genuine divergences in 
perturbation theory. If this was the case, one might expect that the FP ceases to exist as soon 
as one includes in the truncation terms that correspond to nonrenormalizable divergences in 
the perturbative treatment of Einstein's theory. From the point of view of the ERGE, this 
looks quite unlikely, for various reasons. First of all, there exist examples of theories that are 
)erturbatively nonrenormalizable while nonperturbatively renormalizable at a nontrivial FP 
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60] . Second, when one truncates the action Tk, in principle all terms that are retained 



could be equally important. The argument about eliminating terms that contain the Ricci 
tensor only applies when the higher order terms are considered as infinitesimal perturbations 
of the Hilbert action with zero cosmological constant, for then at leading order the on shell 
condition is simply R^ v = 0. However, if the higher order terms are not infinitesimal, the on 
shell condition is much more complicated and there is no indication that they can still be 
eliminated by field redefinitions. In fact for certain classes of terms it is known that they can 



on 



y be eliminated at the price of introducing a number of scalar fields with new interactions 
571 ] . Finally, the calculations reported here already provide evidence to the contrary. In 
the presence of a cosmological term neither R^ V R^ U nor R 2 vanish on shell, so, according to 
the perturbative reasoning, none of the FP's discussed in sections VI and VII should exist. 
If for some reason one is willing to neglect the cosmological constant, in pure gravity the 
first genuine perturbative divergence is cubic in the Riemann tensor, but in the presence of 
matter fields already the one loop logarithmic divergences, which are quadratic in curvature, 
do not vanish on shell. Does this imply that in the presence of matter the FP ceases to 
exist? We have seen explicitly in section VI that at least at one loop this is not the case. 
Admittedly, this is only a partial result, and an "exact" calculation would be necessary to 
definitely settle this point, but it is very strong indication that the Wilsonian approach can 
handle terms that would be troublesome in perturbation theory. For these reasons we also 
believe that nothing special will happen when the Goroff-Sagnotti cubic term (which was 
used to prove that pure gravity is perturbatively nonrenormalizable) will be included in the 
truncation. 

While the systematic derivative expansion cannot for the time being be pushed beyond the 
fourth order, one can still consider different truncations that go beyond the Einstein-Hilbert 
one. So far, the truncation with the greatest number of free parameters that can be dealt 
with exactly is so-called f(R) gravity, where / is a polynomial in the scalar curvature. The 
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calculation of the beta functions of this theory was briefly reported in [ijj] for polynomials 
up to sixth order and has been described in greater detail in section VII for polynomials of 
order up to eight. The most important results of these calculations are the relative stability 
of the results under the increase in the number of terms in the truncation, and the finite 
dimensionality of the critical surface. It appears that the critical exponents do not deviate 
very strongly from the classical dimensions, as expected, so that the terms with six or more 
derivatives are irrelevant. 

Further results involving infinitely many couplings were given in 6jj in the two Killing 
vector reduction of gravity and more recently in 62J in the conformal reduction, where trans- 
verse degrees of freedom of the metric are ignored. There is some work showing that the FP 
exists also in the presence of some types of matter {^j). Some independent evidence in favor 
of a gravitational FP is also found in numerical simulations, both in the causal dynamical 
triangulation approach 64] and in Regge calculus 6^]. This concludes our overview of the 
currently available evidence for a gravitational FP. Let us now make a few comments and 
discuss future prospects. 

The dependence of certain results on the choice of the cutoff scheme (which we called 
"scheme dependence") is sometimes the source of worries. For example, could the FP for 
Newton's constant not disappear if we choose the cutoff function in a perverse way? This 
scheme dependence is the counterpart in the Wilsonian approach of the regularization and 
renormalization scheme ambiguities that are encountered in perturbation theory. As we 
have discussed in section V, such scheme dependence is to be expected for all results that 
concern dimensionful couplings. We have seen that within certain approximations (which we 
called the "perturbative Einstein-Hilbert flow") all the terms with six or more derivatives 
can indeed be made to vanish by a choice of cutoff. However, as we have argued in section V, 
the general properties of the cutoff functions are such that one cannot similarly set to zero 
Newton's constant and the cosmological constant. This is the main reason for the robustness 
of the nontrivial FP. 

Another point that is sometimes a source of misunderstandings is the use of the term 
"exact" in relation to the ERGE, and the "nonperturbative" character of these calculations. 
The beta functions that we have calculated in this paper are rational functions of the cou- 
plings gi. The appearance of the couplings in the denominators suggests that they could be 
regarded as resummations of infinitely many perturbative terms. Thus, the beta functions 
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might still be considered "perturbative" , in the sense that they are analytic in the coupling 
constants. We think that the ERGE is actually capable of obtaining also results that are 
nonanalytic in the coupling constant, but insofar as the FP is present already in the lowest 
order of perturbation theory (see equation (jBJ)) such a degree of sophistication would be 
unnecessary. On the other hand, the beta functions obtained from the ERGE can be said to 
be "nonperturbative" in the sense that their validity is not limited to small couplings. If we 
compare the one loop beta function ([6]) to the beta function (168]) obtained from the ERGE, 
they are seen to have the same form. However, the former calculation by hypothesis is valid 
only for G <C 1; while the latter was obtained from a completely different procedure and is 
not similarly limited. (The difference in the coefficient is not important because the result 
is scheme-dependent anyway.) This suggests that the perturbative result is at least qualita- 
tively valid also for relatively large values of the coupling. Notice that the RG improvement 
modifies the one loop result and produces singularities of the beta functions. The nontrivial 
FP always occurs on the side of the singularity that can be continuously connected to the 
Gaufiian FP. 

The appearance of higher derivative terms in the action at the FP raises the old issue of 
unitarity. The Wilsonian point of view puts that problem in a slightly different perspective. 
A tree-level analysis of the action (128|) plus (1761) shows that it generally contains besides a 
massless spin two graviton other particles with a mass of the order of Planck's mass and 
a negative residue at the pole (ghosts) [58(. Some authors have suggested that, due to 
RG effects, these particles may not propagate 42|, |66|. From a general Wilsonian point of 
view, the presence of a propagator pole at a given mass can only be reliably established by 
considering the effective action at a value of k comparable to that mass. The effective action 
in the FP regime is probably quite different from the effective action at the scale of the 
putative ghost mass, so any conclusion about the mass spectrum based on the FP-effective 
action is probably of little value. This fact is clearly exemplified by QCD. The tree level 
analysis of the QCD action, which is known to be a good description of strong interactions 
at high energies {i.e. near the Gaufiian FP), would suggest a spectrum of particles (quarks 
and gluons) none of which is observed in nature. Since the QCD coupling becomes stronger 
as the energy decreases, the description of strong interactions in terms of quarks and gluons 
becomes extremely complicated long before one reaches the scale of the quark masses, and it 
is believed that if one could actually do such calculations, quarks and gluons would be found 
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not to propagate. It is conceivable that some similar phenomenon may occur in gravity, so 
the FP action (which is the analog of the QCD action) should not be expected to be a good 
guide to the particle spectrum of the theory. The confinement of quarks and gluons is one 
of the outstanding problems in particle physics and it is unfortunately possible that the 
analogous problem in gravity may prove equally hard. This is related to the more general 
question about the low energy action corresponding to the FP action but this lies outside 
the scope of the present paper. We refer the reader to 6^] for some discussions of this point. 

Another remarkable aspect of these calculations is that the one loop flow in the A-G 
plane is essentially the same (aside from nonuniversal numerical coefficients) in the Einstein- 
Hilbert truncation (section IVG) and in the curvature squared truncation (section VI). In 
the latter, A and £ tend logarithmically to zero and the corresponding terms in the action 
diverge. Thus the dynamics becomes dominated by the four derivative terms, while in the 
Einstein-Hilbert truncation it is dominated by the two derivative term. It may be somewhat 
surprising that the structure of the flow should be so similar in spite of such differences in 
the dynamics. This can be at least partly understood by the following argument. In gravity 
at low energies (hence in the perturbative regime) one can consider all couplings to be 
scale-independent, and therefore the relative importance of the terms in the action can be 
determined simply by counting the number of derivatives of the metric. For example, at 
momentum scales p 2 <C Z (recall that Z is the square of the Planck mass), using standard 
arguments of effective field theories, the terms in the action (176]) with four derivatives are 
suppressed relative to the term with two derivatives by a factor p 2 /Z. This is not the case 
in the FP regime. If we consider phenomena occurring at an energy scale p, then also the 
couplings should be evaluated at k « p. But then, if there is a nontrivial FP, Z runs exactly 
as p 2 and therefore both terms are of order p 4 . This argument can be generalized to all 
terms in the derivative expansion (jSJ): at a FP, the running of each coefficient gf- n \ which is 
given by the canonical dimension, exactly matches the number of derivatives in the operator 
0^ n \ so that all terms are of order p A . This may at least in part explain the robustness of 
the results. 

When many terms are taken into account in the truncation, it is hard to have an intu- 
itive feeling for the mechanism that gives rise to the FP. For example, the beta functions 
which are obtained by taking derivatives of (11131) with respect to curvature are exceedingly 
complicated. In fact, they are manipulated by the software and one does not even see them. 
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This is why we have strived in the first few sections of this paper to emphasize the simplest 
approximations. They give a clear and intuitive picture suggesting the emergence of a FP to 
all orders in the derivative expansion. We would therefore like to conclude by overturning a 
common belief: the existence of a nontrivial FP does not require a delicate cancellation of 
terms. The FP appears essentially due to the dimensionful nature of the coupling constants, 
and it can be seen already in the perturbative Einstein-Hilbert flow {i.e. in the approxi- 
mation where one considers just the contribution of gravitons or matter fields with kinetic 
operators of the form —V 2 + E, where E is linear in curvature). More advanced approx- 
imations dress up this simple result with RG improvements and with the contribution of 
additional couplings. The argument in the preceding paragraph suggests that the new cou- 
plings will not qualitatively change the results. And indeed, so far it seems that generically 
such dressing does not spoil the FP. So, to conclude on an optimistic note, one could say that 
it would actually require a special conspiracy by the new terms to undo the perturbative FP. 
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APPENDIX A: TRACE TECHNOLOGY 



The r.h.s. of the ERGE is the trace of a function of a differential operator. To illus- 
trate the methods employed to evaluate such traces, we begin by considering the covariant 
Laplacian in a metric g, — V 2 . If the fields carry a representation of a gauge group G and 
are coupled to gauge fields for G, the covariant derivative V contains also these fields. We 
will denote A = — V 2 1 + E a second order differential operator. E is a linear map acting 
on the spacetime and internal indices of the fields. In our applications to de Sitter space it 
will have the form E = qR 1 where 1 is the identity in the space of the fields and q is a real 
number. 

The trace of a function W of the operator A can be written as 

TyW(A) = W(\i) (Al) 

i 

where A« are the eigenvalues of A. Introducing the Laplace anti-transform W(s) 

POO 

W(z) = / dse- zs W(s) (A2) 



o 



we can rewrite (lAlj) as 

poo 

TrW(A) = / dsTrK(s)W{s) (A3) 
Jo 

where TtK(s) = 'Yl li e~ sXi is the trace of the heat kernel of A. We assume that there are no 
negative and zero eigenvalues; if present, these will have to be dealt with separately. The 
trace of the heat kernel of A has the well-known asymptotic expansion for s — > 0: 

Tr (e"* A ) = — ^ \B (A) s'l + B 2 (A) s'i +l + . . . + B d (A) + B d+2 (A) s + ...1 (A4) 
(47r) 2 l j 

where B n = J d d Xy/gtTh n and b n are linear combinations of curvature tensors and their 
covariant derivatives containing 2n derivatives of the metric. 

Assuming that [A, E] = 0, the heat kernel coefficients of A are related to those of — V 2 

by 

Tre -*(-v2 + E) = _±_ t£t f d d x^trb k (A)E e s k+e ~ 2 . (A5) 
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The first six coefficients have the following form 68]: 

b = 1 (A6) 
R 
~6 



b 2 = ^1-E (A7) 



b4 = lio ( R ^ R ^ ~ RilVR ^ + \ r2 + 6y2jR ) 



+^O^fi^ - ^RE + ^E 2 - ^V 2 E (A8) 

+|e 2 + E 3 + ^ al3 R^ - RTR^ + b -R 2 + 6V 2 i? 
+^O^fi^ + ^EO^O^ + ^EV 2 E - ij„J" 



30 
1 



;^ V2i? + TTc R ^ 2RiiV + ™«? (-642^^^ + A8R^R aP R a f u 



630"' ' 140" A "" ' 7560 v " ~" v " a ~~ ^ 

+ QR.uR^R^ + VTR,J* RJ° RpT ~ ^R\ u p R a p ^R%\)] (A9) 

where fi^ = [V M , V v ] is the curvature of the connection acting on a set of fields in a 
particular representation of the Lorentz and internal gauge group and J M = V a ft a . We 
neglect total derivative terms. The coefficient b§, which is also used in this work, is much 



too long to write here, and can be found in 69]. These coefficients are for unconstrained 



fields. The ones for fields satisfying differential constraints such as h? u and £ M in the field 
decompositions (|46p are given in the following appendix. 

Let us return to equation (1A3|) . If we are interested in the local behavior of the theory 
(i.e. the behavior at length scales much smaller than the typical curvature radius) we can 
use the asymptotic expansion (]A4j) and then evaluate each integral separately. Then we get 

W(A) = — ^- [Qd(W)B (A) + Q i _ 1 (W)B 2 (A) + ... 

( 4vr )2 L > 

+ Q (W)B d (A) + Q^(W)B d+2 (A) + ...], (A10) 

where 

POO 

Qn(W) = / dss- n W(s) . (All) 

JO 

In the case of four dimensional field theories, it is enough to consider integer values of n. 
However, in odd dimensions half-integer values of n are needed and we are also interested 
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in the analytic continuation of results to arbitrary real dimensions. We will therefore need 
expressions for (1A11I) that hold for any real n. 

If we denote the i-th. derivative of W, we have from (IA2D 



POO 

W {i) (z) = (-If / dss^-^Wis) . (A12) 
Jo 

This formula can be extended to the case when i is a real number to define a notion of 
"noninteger derivative" . From this it follows that for any real i 

Qn(W®) = (-lYQn-iiW) . (A13) 



For n a positive integer one can use the definition of the Gamma function to rewrite flAlip 
Mellin transform: 

1 f°° 

Q n (W) = — / dzz n - l W(z) (A14) 

r(n) Jo 

while for m a positive integer or m = 

g_ m (W) = (-l) m iy (m) (0) . (A15) 



More generally, for n a positive real number we can define Q n {W) by equation (1A14|) . while 



for n real and negative we can choose a positive integer k such that n + k > 0; then we can 
write the general formula 

QniW) = T(n + k) j dzzn+k ~ lw(k) ( z ) ■ (A16) 

This reduces to the two cases mentioned above when n is integer. In the case when n is a 
negative half integer n = — 2m + 1 we will set k = m + 1 so that we have 

w 



I 1 \m+l roo 

Q_2_ m±1 (W) = ± ^ dz z~ 1/2 f {m+l \z) (A17) 



Let us now consider some particular integrals that are needed in this paper. As discussed 
in section 1111} there are two natural choices of cutoff function: type I cutoff is a function 
Rk{— V 2 ) such that the modified inverse propagator is Pk{— V 2 ) = —V 2 + Rk{— V 2 ); type II 
cutoff is the same function but its argument is now the entire inverse propagator: Rk(A), 
such that the modified inverse propagator is Pfc(A) = A + Rk(A). 

We now restrict ourselves to the case when E = ql, so that we can write A = —V 2 + ql. 
The evaluation of the r.h.s. of the ERGE reduces to knowledge of the heat kernel coefficients 
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and calculation of integrals of the form Q n y (p^.\y j ■ It is convenient to measure everything 
in units of k 2 . Let us define the dimensionless variable y by z = k 2 y; then Rk{z) = k 2 r(y) 
for some dimensionless function r, Pk(z) = k 2 (y + r(y)) and d t Rk(z) = 2k 2 (r(y) — yr'(y)). 

In general the coefficients Q n (W) wm depend on the details of the cutoff function. How- 
ever, if q = and t = n + 1 they turn out to be independent of the shape of the function. 
Note that they are all dimensionless. For n > 0, as long as r(0) ^ 0: 



Q, 



Similarly, if r(0) ^ and r'(0) is finite, 







o dy 


"i y n 


W 1 ) 


r(n) J 


n (y + r) n 



2 



(A18) 



dtRt 
Pk 



2 . 



(A19) 



Finally, for n 

dtRk 



-m < 



Qn 



pi— m 
r k 



y=0 



dtRk 

pi— m 
r k 



(m) 



£ 

n=0 



[r — yr 



>(«) 



,(m-l) 



y=0 



as [r — yr 



,(n) 



(A20) 

^ r (n+i) _ r (n) _ _^ r (n+i) vanishes at y = 0. This concludes 



the proof that Q n 



pn + 1 



are scheme-independent. 



Regarding the other coefficients Q n 



dtRk 
(Pk+q)' 



we will use the so-called "optimized cutoff function" [3] 

R k ( z ) = (k 2 - z)6(k 2 - z) 



whenever explicit evaluations are necessary, 



(A21) 



With this cutoff dtRk = 2k 2 9(k 2 — z). Since the integrals are all cut off at z — k 2 by the 
theta function in the numerator, we can simply use Pk{z) = k 2 in the integrals. For n > 
we have 



Qn 



dtR 



t-n-fe 



Pk + qYJ n\(l + q) 



_ k 2(n-£+l) 



where q = qk 2 . For n = we have 



Qo 



dtR 



t-n-k 



dtR 



t-n-k 



(Pk + q) e J (Pk + q) 



2=0 

9tR k (z) 



;i + 0) 



-k 2{ ~ l+1) . 



(A22) 



(A23) 



Finally, owing to the fact that the function ^ p *^ + ^y is constant in an open neighborhood 
of z — 0, we have 



Qn 



dtR 



t^k 



(Pk + qY 



for n < 



(A24) 
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This has the remarkable consequence that with the optimized cutoff the trace in the ERGE 
consists of finitely many terms. 
For noninteger n let us calculate 

n (d t R k \ r , -i /2 d t R k { Z ) 



where Pjeiz) = z + (k 2 — z)9(k 2 — z). We change the variable to x — z/k 2 so we have 
fd,R h \ (-i)"+ifc-C-+D f°° <j"+' 2x0(1 - x) 



We find 



so that 



dxx- 1/2 ^-f(x) = 2 

^2 

dxx- l/2 —f{x) = -5 (A27) 



Q-3/2 



5 



We also need some Q-functionals of r^+^p ■ For n > we have 
The function rp^yq^p is equal to ^2^y i n an open neighborhood of z = 0; therefore 



(Pk + qYJ (Pk + qY 



(1 + 9) 

z=0 



;k 2 (- e+ V (A30) 



G-i ( u ) = 7TT^7*~ M , Q n ( /n Rk - f ) = for n < -1 . 



,(p k + q yj (i + qy ' + 

Finally, for the type III cutoff one also needs the following 

Q - {wfw) = ^Stw for " - ; Q " ( wfw) = for " " ■ (A32) 

In conclusion let us address a general problem concerning the choice of the operator O, 
whose eigenfunctions are taken as a basis in the functional space. In some calculations 
the r.h.s. of the ERGE takes the form |Tr 8 ^|^"j^ where A is an operator and q is a 
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constant. Equation (1A10I) tells us how to compute the trace of this function, regarded as a 
function of the operator A + ql. In the derivation of this result it was implicitly assumed 
that O = A + ql. However, the trace must be independent of the choice of basis in the 
functional space. It is instructing to see this explicitly, namely to evaluate the trace regarding 

1 T a t R k (A+ql) 

2 11 P fc (A+gl) 

in general, expanding in q we then have 



^rp r a^fc^i-g±; ag a f unc ti on Q f a. Given any function W(z) we can define W(z) = W(z + q); 



1 f°° 

Q n {W) = — - / dzz n - l W{z + q) 
1 \ n ) Jo 

= — - / dzz n - 1 (W(z) + qW'(z) + -q 2 W"(z) + -q 3 W"'(z)...) 
T(n) J 2! 3! 



QniW) + qQ n {W') + -q 2 Qn(W- 



3! 



q 3 Qn(W") + ... 



Qn(W) - qQn-xiW) + Ig 2 Q n _ 2 (W0 - ^Q n _,(W) . . . 



(A33) 



where in the last step we have used equation (lA13h . Using (1A10I) for the function W we 
then have 

W[A] 



(47T 
1 



(47T 



Qa (W)Bo(A) + Qi_ 1 (W)B 2 (A) + ... + Q (W)B 2d (A) - 
(q%(W) - qQi^W) + Ig 2 Q f _ 2 (W/) - Ig 3 g f _ 3 (H/) 



<?-_! (W) - qQ!_ 2 (W) + lg 2 Q f _ 3 (H/) - lg 3 g f _ 4 (H/) 



B (A) 
B 2 (A) 



Qo(W) - qQ-^W) + ^q 2 Q~ 2 (W) - ^q 3 Q^(W) + ...) B 2d (A) 



+ 



(A34) 



We can now collect the terms that have the same Q-functions. They correspond to the anti- 
diagonal lines in flA34j) . Using equation (IA5I) one recognizes that the coefficient of Q±_ k is 
B 2 k(A + ql). Therefore 

1 



TrW[A] 



(47T) 



2 L 



Q*(W)B Q (A + ql) + Q i+1 (W)B 2 (A + ql) 



(A35) 



+ ... + Q (W)B 2d (A + ql) + .. 

which coincides term by term with the expansion of TrVU[A + g] using the basis of eigenfunc- 
tions of the operator O = A + ql. This provides an explicit check, at least in this particular 
example, that the trace of this function is independent of the basis in the functional space. 
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APPENDIX B: SPECTRAL GEOMETRY OF DIFFERENTIALLY CON- 
STRAINED FIELDS 



In this appendix we work on a sphere. Consider the decomposition of a vector field 
into its transverse and longitudinal parts: 

A, - Al + 



The spectrum of —V 2 on vectors is the disjoint union of the spectrum on transverse and 

R 
d 



longitudinal vectors. The latter can be related to the spectrum of —V 2 — 4 on scalars using 



the formula 

- V 2 V,$ = -V, (V + $. (Bl) 
Therefore one can write for the heat kernel 

Tre-<- v2 ) \ A = Tre-<- v2 ) \ A r +^e~<-^) |* (B2) 

The last term has to be subtracted because a constant scalar is an eigenfunction of —V 2 — 4 
with negative eigenvalue, but does not correspond to an eigenfunction of —V 2 on vectors. 
The spectrum of —V 2 on scalars and transverse vectors is obtained from the representation 
theory of SO(d + 1) and is reported in table IVIIl 

A similar argument works for symmetric tensors, when using the decomposition (1461) . 
One can use equation 

and equation 

- V 2 (v„V„ - - d 9^ 2 ^j ° = (v^V, - - d 9,uV 2 ^j (-V 2 - j^jR^ ° (B4) 

to relate the spectrum of various operators on vectors and scalars to the spectrum of —V 2 
on tensors. One has to observe that the d(d + l)/2 Killing vectors are eigenvectors of 
— V 2 — jjdZYjR on vectors but give a vanishing tensor h^, so they do not contribute to the 
spectrum of —V 2 on tensors. Likewise, a constant scalar and the d + 1 scalars proportional 
to the Cartesian coordinates of the embedding R n , which correspond to two the lowest 
eigenvalues of —V 2 — -^jR, also do not contribute to the spectrum of tensors. So one has 
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for the heat kernel on tensors 




The last exponentials can be expanded in Taylor series as Y2m=o c m R m and these terms 
can be viewed as modifications of the heat kernel coefficients of —V 2 acting on the differen- 
tially constrained fields. To see where these modifications enter, recall that the volume of 
the sphere is 

so that 

J d d x^gtih n oc TT^ . (B7) 

This means a coefficent c m from the Taylor series will contribute to a heat kernel coefficient 
for which 2m = n — d. So there are contributions to b n only for n > d. 

We have discussed how the negative and zero modes from constrained scalar and vector 
fields affect the heat kernel coefficients of the decomposed vector and tensor fields. These 
modes have to be excluded also from the traces over the constrained fields; this is denoted 
by one or two primes, depending on the number of excluded modes. This can be done 
by calculating the trace and subtracting the contributions to the operator trace from the 
excluded modes. Thus the trace with m primes is 

m 

Tr'-' [W(-V 2 )} = Tr [^(-V 2 )] - £ D l (d, s) W (A, (d, s)) (B8) 

1=1 

where Xi(d,s) are the eigenvalues, Di(d,s) their multiplicities, both depending on the di- 
mension d and on the spin of the field, s. The eigenvalues and multiplicities for the m-th 
mode of the Laplacian on the sphere are given in table IVIIi 

The expressions that we will need are those for the cases where one mode is excluded 
from the transverse vector trace (s = 1, m — 1), or one or two modes from the scalar trace 
(s = 0, m — 1, 2), each one in two and four dimensions. The results obtained by calculating 
the corresponding multiplicity and eigenvalue from table IVIII are given in table IVIIII To see 
what is the relevant contribution to one of the heat-kernel coefficients, one can expand the 
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Spin s 



Eigenvalue Xi(d, s) 



Multiplicity Di(d,s) 



l(l+d-l) p. , _ n i 
d(d-l) l 



(2i+ri-l)(i+d-2)! 



d(rf-l) -ft, t — 1, Z , 



f(Z+rf— l)(2Z+rf— 3)! 
(d-2)!(Z+l)! 



*(Z+rf-l)-2 p. , _ 9 o 
rf(rf-l) 1 ~ z ' ° ' 



-l)(rf— 2)(Z+rf)(Z— l)(2Z+rf— l)(Z+rf— 3)! 
2(d-l)!(*+l)! 



TABLE VII: Eigenvalues and their multiplicities of the Laplacian on the d-sphere 





8=1 


s=0 


m = 1, d = 2 


3W(f) 


W(0) 


m = 1, d = 4 


10W(|) 


W(0) 


7Ti = 2, d = 2 




W(0) +3W( J R) 


7Ti = 2, d = 4 




W(0) + 5W(|) 



TABLE VIII: £™i A (d, s) W (Aj (d, a)) for s = 0, 1, d = 2, 4, m = 1, 2 
obtained expression in R. For the case s = 0, d = 4, m = 2 one has for example 

2 



^ A (4, 0) W (A/ (4, 0)) = (0) + 5iy (j 
J dz^ (w (0) + ^W" (0) 



' R 2 W" (0) + — ^— R 3 W" (0) + . . . )(B9) 



= 4Wi ClXv/ "r iUj + ^ H/ ' lUj+ 108^ ' ~ ' 36-27- 
From this one sees that, in this case, the b 2n receive a correction for n > 2. In two 
dimensions, that would be already the case for n > 1. The full list of heat kernel coefficients 
of —V 2 in 4d is given in table IIXI 



APPENDIX C: PROPER TIME ERGE 



Let us start from the ERGE for gravity in the Einstein-Hilbert truncation with a type 
III cutoff, written in equation (ISTJj) . Define the functions: 

d t R k (z) _ . . R k (z) 



A k (z) 



BJz) 



CM = . (ci) 



z + R k (z) " iV "' z + P k (z) z + R k (z) ' 

The term in equation (1601) containing C is nontrivial. To rewrite it in a managable form we 
take the Laplace transform: 



C k (z)= I dsC k (s)e 
lo 



(C2) 
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V 

V 


VT 


T 


TT 


TTT 

_L ± _L 


trb 


1 


A 


Q 
O 


1 n 

1U 


Q 


c; 



11 L>2 


R 
6 


2R 
3 


H 
4 


5R 
3 


3R 
2 


5R 
6 


trYtA 


29R 2 
2160 


43iJ 2 
1080 


7R 2 
1440 


11_R 2 
216 


81R 2 
2160 


-R 2 
432 


trb 6 


37R 3 


R 3 


541R 3 


1343.R 3 


-319i? 3 


311R 3 


54432 


17010 


362880 


136080 


30240 


54432 


trb 8 


149R 4 


2039R 4 


157i? 4 


2999.R 4 


683R 4 


109R 4 


6531840 


13063680 


2488320 


3265920 


725760 


1306368 



TABLE IX: Heat kernel coefficients for S 4 . The columns for the transverse vector (VT) and 
transverse traceless tensor (TTT) are obtained from equations (lB2j) and (|B5|) in d = 4. Note that 
the excluded modes contribute to trb n only for n > 4. 



Since the operator dt{^2 — 2 A) commutes with A2 — 2 A, we can write 



C fe (A 2 -2A)d t (A 2 -2A) 



dsC k (s) d t (A 2 -2A)e 



-s(A 2 -2A) 



°°-C k (s)d t e- s ^- 2A K 
s 



(C3) 



Laplace transforming also A k and B k , the first term in equation fl60l) becomes 



ds 



A k (s) + B k (s)r 1 --C k (s)d t 
s 



Tre -s(A 2 -2A) _ 



(C4) 



This is the functional RG equation in "proper time" form {71, Q, H]- Note that the first 
term corresponds precisely to the one loop approximation. The trace of the heat kernel can 
be expanded 

1 



Tre -s(A 2 -2A) _ e -s(-2A) 

whereas for the ghosts 



( 47r )d/2 
1 

(4^2 



dxsfgti 
dx^fg 



d 




IS 2 + 









+1 



0(i? 2 



d d(7 — 5d)^ d,-. 



12 



O^ 2 



(47r)rf/2 
1 

(47r) d / 2 
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SO 



The ERGE then takes the form: 



t± k 



( 47r )d/2 

"d(7-5d) 



(d(d+l) 



4 



(A k + r/fifc - 2atAG fc ) - dQ ± (A k 



Qi-x (A k + rjB k -2d t AC k ) 



d + 6 



Q±^ (AO 



R + 0(R 2 



where W(z) = W(z — 2A). Using an optimized cutoff one can now reproduce equations ( |62|) 
and fl63|) . However, in this way the sums in equation (16TT) can be resummed for any type of 
cutoff shape. 



APPENDIX D: CUTOFF OF TYPE IB WITHOUT FIELD REDEFINITIONS 



We collect here the formulae for the beta functions of A and G in the Einstein-Hilbert 
truncation, using a cutoff of type lb and without redefining the fields £ M and a. The ERGE, 
including the contributions of the Jacobians, is 

dT k 
dt 



1 T d t Rk + yRk 
2 r(2) R,_2A + 4#^R 



d(d-i) 

1 d t R k (2P k + ^P - 2A) + V R k (P k + z + ^R - 2A) - 2d t AR k 
+ 2 r « (P fc -f)(P fc + ^i?-2A) 

<9 t P fe + T]R k 1 , 



+-Tr 



2 ir(0) P fc - 2 A + '±?R + 2 Tr W (P fc + ±£R - 2 A) P fc (P fe - 

P /d-4 



x 



dtRk 
+77 (P fc 3 -z 3 ) + 



# d 2 - 6d + 4 
3P, 2 + 2P fc ( ~ ; — ^-P - 2A 



-2&a 



d(d- r 

d 2 - 6d + 4 
d(d- 1) 

Pi 



P - 2A ) (P; 



P /d-4 



d- 1 



-Tr (1) 



d- 1 V d 
d t R k Tf/ 2(P fc -f)9 t P fc 



P - 2A P 



p H 



1 



dtRk 



-Tr'/, 



2 ' p fc - 2(£iy) ^-Rfc 



(0)- 



Tr 



which gives 

d£k 
dt 



dxy/g{ ( 

+vQ. 

-d t AQ 



P k(Pk ~ sri) 

(d+l) ((d + 2) P fc -4A+) d t P fc 
4P fc (P fc — 2 A) 

d (d + 1) P fc 2 + 2d (2 - 2A) Pfc + 2z (z - 2A)) R k 
4Pfc 2 (Pfc -2 A) 

(dPfc + z) Pfc 



(0) (ft-f)^ 
(0) ~P~ 



(Dl) 



Pfc 2 (Pfc -2 A) 
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+R 



-Q, 



[(d 4 - 2d 3 - Id 2 + 18d - 16)P 2 + 8(d 2 - d + l)A(P fc - A)] d t R k 



4d(d-l)P 2 (P fc -2A) 2 
+ V Q* ( (d 4 - 4d 3 + 3d 2 + 12d - 16)P 3 + 2(d - 1) 2 (2A + (d - 4)z)P 2 
+ (2(d 2 - 6d + A)z 2 + 4(d 2 -d + l)zA + 8(d - 1) 2 A 2 ) P fc 
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+d t A Qa 



+Q 
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-d t A Q, 
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With the optimized cutoff one obtains the coefficients 

d 4 + 5d 3 + 4d 2 + 4d + 16 + (4d 3 + 12d 2 - 32d - 32)A 



(D2) 



A 1 
A 2 

A 3 
Pi 



-2tt- 



(47r)lr(f + 3)(2A-1) 



16tt 



d 3 + 13d 2 + 48d + 28 - (4d 2 + 44d + 120)A 
(4vr)l(d + 2)(d + 4)(d + 6)r(|)(l - 2A) 
647r(d + 5) 



(47r)2(d 2 + 6d + 8)r(f)(l-2A) 



Bo 
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and the beta functions in four dimensions 
H = -2A 



1 (21 - 105A + 142A 2 - 78A 3 + 92A 4 )G + 3334-13309A+12416A 2 -780A3 ^ 2 
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